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ABSTRACT 


A 2-parameter forecasting model which includes the effects of vertical vorticity advection and turning of the 
vortex tubes is briefly described. Contributions of the above-mentioned terms are discussed and an example is 


presented. 


The importance of consistent truncation over a large grid is pointed out. 


1. INTRODUCTION 


In the summer of 1957, the Joint Numerical Weather 
Prediction Unit (JNWP) changed from an IBM 791 to 
in IBM 704 electronic computer. Because of the higher 
capacity of the latter, the forecast area was extended to 
cover the Northern Hemisphere north of latitude 13°. 
Rather than recode the two-level, geostrophic model in use 
prior to the change of machine, it was decided to test a 
‘lightly more elaborate two-level model in the hope of im- 
proving the forecasts. This model was designed by Lt. Col. 
Philip D. Thompson, then Chief of the Research and De- 
lopment Section of JNWP [1]. 

In testing this model some difficulties were encountered 
which had not been fully anticipated. These difficulties 
vere due in part to attempted simple extension of the model 
from the one previously used, but more importantly to 
the inereased importance of accumulated systematic small 
rors over a large grid. 

This report is concerned mainly with consequences re- 
sulting from including the vertical advection of vorticity 


id the twisting terms in the prognostic equations. 


x, Present affiliation: Navy Numerical Weather Problems Group, U.S. 
“aval Postgraduate School, Monterey, Calif. 


2. DESCRIPTION OF THE MODEL 
The model’s information levels are at 400 and 800 mb. 


As is customary, the -profile (4) is assumed to 


have a smooth shape, and is for this particular model 
symmetrical in the vertical about p=600 mb. and zero for 
p=200 and p= 1000 mb. 


Application of the vorticity equation to the information 
levels results in two predictive equations: 
On, Ov; 


One Ov, 


where the subscripts 1 and 2 refer to 400 and 800 mb., 
respectively; notations are conventional. 
We define the new variables 


n=) (m—n2); (vi—v2) 
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Addition and subtraction of (1) and (2), utilizing (3), 


result in 
on, - ( Ow 


Om_ 1 ( Ov; = 
+5 | op 2 op +5 k Vu, x op Va. >< op =(. 


These equations are simplified when we observe that 
Where the proportionality factor is 
assumed to be less than unity, and that the model can not 
distinguish between the p-derivatives of v and 7 at the 
two information levels. Asa further step toward simplifi- 
cation, the derivatives with respect to p are replaced by 
finite differences as follows : 


Ow, 
op op P op 


on _ On, 2n!' ov; 
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Where P=400 mb., and w refers to the 600-mb. surface. 
The prediction equations for 7 and 7' now take the form: 


2Aw 


Vo Xy'= 0 (5) 


It is pointed out that the fourth term in (5), resulting 
from the divergence terms of (1) and (2), is of the 
same form and magnitude as the fifth term in (5), which 
is the contribution of the vertical advection of vorticity. 
This implies that if the vertical advection and twisting 
terms were omitted from the vorticity equation, one 
should approximate the divergence term in such a way 
that similar terms are not introduced. As will appear 
later, this is of some importance for the vorticity balance. 
It follows from the symmetry imposed upon the o-dis- 


tribution that the wind vector v is divergence-free and can 
therefore be derived from a stream function y. This is 
accomplished by solving the balance equation. 


In deriving the differential wind and vorticity v' and 7’, 
respectively, the geostrophic approximation is invoked, 
leading to the expressions: 


v=LExvh; (7) 


2f 


where A is the thickness of the layer 800-400 mb. Since 
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v=kxvy and 
the two prediction equations (5) and (6) contain the 
three unknowns y, h, and w. The third equation needed 
for making the system of equations complete is the 


adiabatic equation, applied at 600 mb., 


oh o 
where 
oo 


and p is the density and 0 the potential temperature. 


Utilizing the relations (7), (8), and (9), we may write 
the predictive equations (5) and (6) as follows: 


_1+2A gw gA 


2 Vo'Vh=0 (11) 


In solving the complete system (9), (11), and (12), one 
proceeds as follows : 


1. solve equation (12) for O//d¢ and extrapolate in 


time; 
2. compute » from (9), using 0A/d¢ from step 1; 


3. solve (11) for y at time ¢+ Af; 
4. repeat the steps 1 to 3. 


It can be disputed whether it is worth while including 
small terms like the twisting and vertical advection terms 
in such a crude model. As already mentioned, the model 
cannot distinguish between the vertical advections at the 
two information levels and the same applies to the twist- 
ing term. 

In testing the model, however, it was decided to study 
the effects of the various terms in equation (11), and as 
a result three versions of (11) were dealt with: these are: 


1. Equation (11) is replaced by 
13) 
ot — {P V*h=0 (13) 


i.e., the third and fifth terms in (11) have been omitted 
and A has been set equal to 4; 
2. Equation (11) is replaced by 


i.e., the twisting has been omitted ; 


int 


is 


| 
iJ 
| te 
to 
Op th 
al 
In 
he 
Of 
th 
sy 
ar 
by 
diy 
cu 
i, 
| pel 


‘ite 


11) 


12) 


me 


13) 


14) 


Apart 1959 


3, Equation (11) is replaced by 


_ ga VA— 


iP -Va.Vh=0 


(15) 


ie, the term . V°A, originating from the divergence 


ierms in (1) and (2), has been neglected, and A set equal 
to unity. The latter is inconsistent with the assumption 
that A be less than unity, but this inconsistency is not of 
any significance. 


3. COMMENTS ON VORTICITY BALANCE 
AND FICTITIOUS SOURCES OF VORTICITY 


Omitting friction, we may write the vorticity equation 
as follows: 
Ov 


divn Xe —=0. 


op (16) 


Integration of (16) over an area A of a pressure surface 
hounded by a curve C gives 


1 OV 4 


where d.1 is an area element, n a unit vector normal to 
the unit vector t which in turn is tangential to the curve C 
of which dC is an element; t, n, and k form a right-hand 
system. 


We substitute for vy from 


v=kXV¥+VX (18) 
into (17) which takes the form 


Thus the mean vorticity ¢ may be changed as follows: 

(a) by advection of vorticity across the lateral bound- 
wies by the divergence-free part of the wind; 

(b) by vorticity advection across the lateral boundaries 
by the non-rotational part of the wind and due to its 
livergence (VX. and ydivv) ; 

(¢) vertical advection of circulation along the boundary 
curve 

If we impose the boundary condition v-n=0 and o=0, 
«, if the system is closed, the mean vorticity is inde- 
pendent of time. Hence, as ons as the vorticity equation 


is not approximated there are no vorticity sources. This 
‘Sof interest to bear in mind when dealing with simple 
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models where approximations are inevitable. For ex- 
ample, in omitting the advection of vorticity by the non- 
rotational part of the wind in the vorticity equation and 
keeping the divergence term, a false vorticity source may 
thereby be introduced, upsetting the vorticity balance. 
This argument applies also to the twisting and vertical 
advection terms in the vorticity equation; retaining one 
of them and omitting the other introduces a spurious 
vorticity source. 

As will be shown later, this introduction of spurious 
vorticity sources is not to be taken lightly when dealing 
with forecasts over a large area. For a hemispheric grid 
it introduces a large-scale error devastating for the use- 
fulness of the forecast. 


4. DISCUSSION OF THE FORECASTS 


All forecasts discussed here were made using the 
JNWP octagonal grid which covers the Northern Hemi- 
sphere north of latitude 13°. 

In investigating the performance of the models, the 
contributions of certain of the terms of the prediction 
equation for were computed separately and then 
accumulated for 12-hour periods. 

1. Forecasts using equation (13) were carried to 48 
hours for 4 cases with results which were similar with 
respect to the features of interest here. In these four 
cases the contribution of the third term of (13) was 
isolated. The forecast with initial time 1500 er, April 
3, 1958, may be taken as representative of this group, and 
is reproduced as figure 1. During the first 12 ours the 
vertical advection of vorticity contributed a negative 
height change everywhere with maximum values of 
—600 feet. During the second 12-hour period an area 
of positive changes up to 160 feet appeared, but decreased 
negative changes persisted over most of the map. During 
the third 12-hour period some small negative changes 
persisted, but the positive contributions greatly enlarged 
with maximum values up to 370 feet. During the last 
12-hour period the contributions were everywhere positive 
with maximum values of 400 feet. The contribution of 
this term appeared to be mainly large-scale with little 
detail on the scale of the principal synoptic features. 


The phase relationships between A and y and between 
A and w are illustrated in figure 2. The vertical velocity 
field, initially as well as throughout the 48-hour forecast 
period, was well ahead of the streamline field (upward 
motion ahead of a trough). On the other hand, the thick- 
ness field started out lagging behind the stream function, 
caught up at around 24 hours, and was well ahead at 48 


2The quantity y occurring in the figures has the dimension length and 
is derived from ¥ through the relation y=“ yw where g is the accelera- 


fe 
tion of gravity and fe is the Coriolis parameter at 45° latitude. In the 
captions to the figures, y is alternatively referred to as height and stream 


function. 
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Figure 1.—Change in 500-mb. height due to vertical advection of vorticity, labeled in tens of feet. 


hours. The error maps show that the stream function vorticity as time increased. 
field was moved too slowly, and indicate too fast a move- 2. A forecast using (14) was carried to 48 hours, aga! @ |S. 
ment of the thickness field. The w-field versus A-field using as initial data 1500 amr, April 3, 1958. tri 
shows that w and h were positively correlated initially In this forecast the thickness field kept its lag better th 
and negatively correlated after 24 hours. This would than with the previous model, but during the 36—48-how les 
seem to explain the reversal in the vertical advection of period it caught up with the stream function field" * 
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Ficure 2.—-Phase lags between thickness and stream function and thickness and vertical velocity. 
lines indicate ridges and short solid lines troughs in the h-field (not contoured); the labeled solid lines are y-contours. 


In portraying the h, ¥ phase, dashed 
In the two 


maps showing the phase between thickness and vertical velocity (h,w), solid lines are thickness contours (in tens of feet) and dashed 


lines contours for vertical velocity in em. see.~!. 


parts of the map. The vertical advection contribution 
8 shown in figure 3. There was no reversal in the con- 
tribution of the vertical advection of vorticity term, but 
the negative contribution to stream function change was 
ess during the 24-48-hour periods than during the first 
“hours. Vertical velocity had much the same phase 


lationship to stream function as before, but there is con- 


siderable difference in detail between the fields. The 
24-hour vertical velocity with (14) was more intense, and 
the 48-hour vertical velocity slightly more intense than 
results from using (13). 

3. Six forecasts were carried out using equation (15) 
as one of the basic prognostic equations. Here again the 
results were similar with respect to the features under 
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Figure 3.—Change in 500-mb. height due to vertical advection of vorticity, labeled in tens of feet. 


consideration, and the case of 1500 emt, April 3, 1958, 
may be taken as typical. In this case the last two terms 
of equation (15) were accumulated separately for each 
12 hours, the results being given in figure 4. It is seen 
that the integrated values almost everywhere have oppo- 


site signs and tended to cancel each other. The vertical 
advection term was the larger of the two but this } 
mainly because of different degrees of truncation. I 
the vertical advection term, the Laplacian of thicknes 
is formed using the mesh length of d, whereas both of the 
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Figure 4.—Change in 500-mb. height due to vertical advection (solid lines) and twisting (dashed lines) terms. 


first derivatives in the twisting terms have 2d as mesh 
length. When the inconsistency of truncation is removed 
by combining the two terms it becomes clear, as pointed 
out earlier, that the net contribution over the grid of the 
‘wo terms depends only on conditions at the lateral 
boundaries. 

When the last two terms of (15) are combined before 
finite differences are taken, the accumulated sums take on 
the Values given in figure 5. These may be compared with 
the results shown in figure 6, which is the sum of the 
‘wisting and vertical advection terms of figure 4. These 


results should be the same except for truncation errors. 
It thus appears that on a large grid inconsistent truncation 
of small but systematic terms may alone account for sig- 
nificant forecast error. 

It is also to be remarked that in this model and informa- 
tion grid the vertical velocity terms contributed very little 
in comparison with the error remaining. However, it 
would be unwise to conclude from this experiment alone 
that the terms are really unimportant either to the atmos- 
phere or to more sophisticated models. 

This model, when corrected for the systematic error 
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Figure 5.—Change in 500-mb. height due to combined vertical advection and twisting terms consistently truncated. 


described above, showed a tendency to move the sub- 5. CONCLUSIONS ti 
tropical Highs erroneously westward, and a tendency to- The conclusions reached in this study may be suit 

ward overdevelopment of Lows. The spurious westward — jyarized as follows: t 
movement is presumably a manifestation of this model’s 1. Each of the vertical advection and twisting term ‘ 
incapability of properly forecasting the long waves as jn the vorticity equation tends to have a uniform sigh , 


discussed by Wolff [2] and Cressman [3]. The over-  oyer the entire grid. In a closed system their net col 
development is attributed to the third term in equations tribution is, however, zero when taken together. Omit- 
(14) and (15). This term would seem to be incapable ting one and retaining the other of these two terms leads 
of properly accounting for observed genesis of mean to spurious creation or destruction of vorticity, devas 
vorticity ¢. tating to a hemispheric forecast. 
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Figure 6.—Change in 500-mb. height due to vertical advection and twisting terms inconsistently truncated. 


2, When both terms are included, inconsistent trunca- REFERENCES 
tion may lead to considerable error in the forecast. 1. P. D. Thompson, “A Two-Level Model with Effects of Vertical 
m- 3. In approximating the divergence term in the vor- Velocity Advection, Irregular Terrain, and Variable Static 


Stability,” Office Note No. 6, Joint Numerical Weather Pre- 


ticity equation, precaution should be taken not to intro- diction Unit. June 1987 


. 
‘ duce a small but systematic error of the type referred to 9 Pp. M. Wolff, “The Error in Numerical Forecasts Due to Ret- 
4 under conclusion 1. rogression of Ultra-Long Waves,” Monthly Weather Review, 
it 4. The contribution of the vertical advection and twist- vol. 86, No. 7, July 1958, pp. 245-250. 

in win 3. George P. Cressman, Barotropic Divergence and Very Long 
d g ; iS when taken oget en nt m Atmospheric Waves,” Monthly Weather Review, vol. 86, No. 
as- model in the six cases tested. 8, August 1958, pp. 293-297. 
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_ New Weather Bureau Publications 


Technical Paper Series 


Numbers issued in the Technical Paper Series since 
the last publication of the list (MWR Aug. 1957) are 
/ given below. All are for sale by the Superintendent of 
Documents, Government Printing Office, Washington 
25, D.C. 


No. 


No. 


No. 


No. 


15. “Maximum Station Precipitation for 1, 2, 3, 6, 12, 
and 24 Hours”: Part X XI, Illinois, 1958, 50 cents; 
Part XXII, Ohio, 1958, 65 cents; Part X XIII, Cali- 
fornia, 1959, $1.50. 

29. “Rainfall Intensity-Frequency Regime”: Part 2, 
Southeastern United States, 1958, $1.25; Part 3, The 
Middle Atlantic Region, 1958, 30 cents. 

32. “Upper-Air Climatology of the United States”: 
Part 1—Averages for Isobaric Surfaces, Height, 
Temperature, Humidity, and Density, 1957, $1.25; 
Part 2—Extremes and Standard Deviations of Aver- 
age Heights and Temperatures, 1958, 65 cents; Part 
3—Vector Winds and Shear, 1959, 50 cents. 

Part 1 contains tabulations and charts by months 
covering the Caribbean area and the United States 
including Alaska for the period 1946-55 and for the 
surface to the highest level for which sufficient data 
were available. Part 2 contains the monthly ex- 
tremes and standard deviations for the temperature 
and height data presented in part 1. Part 3 presents 
seasonal charts of mean vector winds and standard 
deviations covering only the contiguous 48 States and 
for various pressure surfaces from 700 mb. to 30 mb. 
33. “Rainfall and Floods of April, May, and June, 
1957 in the South-Central States,” 1958, $1.75. 


Contains analyses of the meteorological situation 
from the general circulation, synoptic, and radar 
standpoints, a description of the flooding and dam. 
age, and tables of the daily and hourly rainfall at 
recording and nonrecording stations in the area. 

No. 34. “Upper Wind Distribution Statistical Parameter 
Estimates,” 1958, 40 cents. 

Tabulations of such parameters of the wind data 
as direction and magnitude of vector mean wind, 
standard vector deviation, mean and standard devia- 
tion of zonal and meridional components, correlation 
of meridional to zonal components, angle of rotation 
of elliptic axes of wind distribution, standard deyvia- 
tion of components along major and minor axes of 
ellipse, ratio of standard vector deviation to square 
root of difference between squares of scalar and 
vector means. Covers Northern Hemisphere and 
mostly period 1948 to 1953. 

35. “Climatology and Weather Services of the St. 
Lawrence Seaway and Great Lakes,” 1959, 45 cents. 

Descriptive climatology of area with illustrative 

charts. 


No. 


Storm Data 

A new publication Storm Data was begun in January 
1959. This monthly periodical consists of the informa- 
tion on storms and storm damage formerly contained in 
Climatological Data—National Summary. Tt has been 
separated from the parent publication and the tabula- 
tion somewhat streamlined to speed up presentation of 
the storm statistics. Price: 15 cents per copy: annual 


subscription $1.50 domestic, foreign. 
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SOUNDINGS OF TERRESTRIAL RADIATION FLUX OVER WISCONSIN’ 


P. M. KUHN, Office of Meteorological Research, U.S. Weather Bureau 
V. E. SUOMI and G. L. DARKOW, Department of Meteorology, The University of Wisconsin 
(Manuscript received January 6, 1959; revised April 7, 1959] 
ABSTRACT 


Results and conclusions drawn from 50 radiometer-sonde ascents over Wisconsin are presented. Atmospheric 


infrared radiation flux is separated into upward, downward, and net components. 


Radiation divergence and 


resulting cooling are computed. Seasonal averages of these direct observations are compared. 


1. INTRODUCTION 


In the past, direct measurements of the atmosphere’s 
radiative balance have been costly and complex. Suomi, 
Staley, and Kuhn [1] recently reported on a radiometer- 
snde which was designed as a routine technique for ob- 
taining vertical profiles of the flux of infrared radiation 
inthe atmosphere. This instrument has been successfully 
adapted to balloon-borne flights, and briefly, consists of 
the economical net radiometer described by Suomi and 
Kuhn [2] attached to the standard U.S. Weather Bureau 
radiosonde. An electrically driven sequencer and simple 
cireuit modifications are necessary. 


SYMBOLS USED IN TEXT 


o7/0¢ Change in temperature in °C. per day. 

y Acceleration due to gravity, ¢.g.s. units. 

Rye Net radiation in langleys per minute at top of 
layer (/2yy>0 directed upward). 

Rvp Net radiation in langleys per minute at bottom 
of layer (Ry_>0 directed upward). 

Cp Specific heat of air, c.g.s. units. 

Ap Thickness of layer in millibars. 

Stefan-Boltzmann constant, ¢.g.s. units. 

ry Observed temperature of bottom sensor. 

I’, Observed temperature of top sensor. 

( Internal conduction for radiometer in langleys 
per minute. 

(, Top conduction for radiometer in langleys per 
minute. 

(, Bottom conduction for radiometer in langleys per 
minute. 

Ey Error term for bottom half of radiometer, lang- 
leys per minute. 

E Error term for top half of radiometer, langleys 
per minute, 

Nore. Values of radiation flux or net radiation >0 indicate a 


fux or net flux directed away from the air-ground interface. 


‘The research reported in this study is supported by the Office of Me- 
‘orological Research, U.S. Weather Bureau. 


2. TYPICAL RADIOMETER FLIGHTS 


The radiometer-sonde flights were carried out at the Uni- 
versity of Wisconsin at Madison, Wis. The flight program 
included (i) direct measurement of the upward, down- 
ward, net, and total components of infrared radiation flux ; 
and (ii) computation of radiative cooling of the 
atmosphere. 

Figures 1 through 3 show three typical nocturnal 
ascents. In each, the upward, RA, downward, RY, and 
net, Ay, radiation currents in langleys per minute are 
shown, together with the air and dewpoint temperatures. 
The radiative temperature change of the air column is 
calculated from the vertical divergence of radiation in 50- 
mb. layers up to the 200-mb. level and for 25-mb. layers 
above that from the expression, 


oT  —g(Ryv—Ryx) (1) 


ot c,Ap 


This profile, given in degrees Celsius change per day, 
appears on the right for each figure. 

Figure 1 is a clear summer night sounding and illus- 
trates characteristic features of many similar flights. In 
all such cases relative warming occurs just beneath any 
inversion, including the tropopause. In particular, the 
warming was evident at 815 mb. during the flight of 
July 7, 1958, beneath a low-level inversion. Again at 
270 mb., warming occurred beneath the tropopause in- 
version. The latter phenomenon is clearly visible im- 
mediately below the tropopause in the summer mean 
cooling profile. Frequently, such radiative warming be- 
comes positive, at times reaching from 2° to 3° Celsius 
heating per day. For example, a warming of 1.0° C. per 
day occurred at 270 mb. on the July 7, 1958, ascent. 

Secondly, in many instances we find the strongest cool- 
ing just above inversions or cloud layers. In detail, this 
may be seen in figure 1 at 680 and 165 mb., reaching over 
2.0° C. in each case. Figure 1 might be compared with 
the observed cooling rate obtained by Brewer and Hough- 
ton [3] for July 20 and 21, 1956. 
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Figure 1.—Clear summer radiometer-sonde ascent. 


The clear nocturnal winter sounding, figure 2, pro- 
vides a typical cooling profile often observed in the winter. 
The average cooling of this winter profile, except for the 
microlayer, is 1.0° C. per day, as compared with a sum- 
mer average cooling rate of 1.7° C. per day. Thus direct 
radiation measurements allow a computation of radiative 
cooling far above the limits of chart evaluation. Between 
250 and 200 mb., there is indication of the ascent passing 
above effective atmospheric radiators. Both the upward 
and downward radiation currents approach a constant 
value even though there is a positive gradient in the tem- 
perature profile. The last plotted point in the upward 
stream is doubtful. 

A cloudy winter ascent, figure 3, demonstrates the radia- 
tive effects of cloud. Below the stratocumulus base at 
870 mb. there is radiative warming. At the cloud top, 
however, there is a sharp gradient in the downward 


00 10 20 30 40 50 -70° -60° -50° -40° -30° -20° -0° Oo 20° -5°-2" 2s 
TEMPERATURE (‘C) 


° 
Yay 
July 7, 1958, 2000 cst, Madison, Wis. 


directed radiation, RY, resulting in strong radiative d- 
vergence and a cooling rate of 8.5° C. per day. The 
break into the clear night sky at 725 mb. sharply reduees 
the downward propagating current. The gradient of RY 
increases from —0.03 ly./min. per 50 mb. to —0.09 ly, 
min. per 50 mb. In the layer between 100 and 150 mb, 
regardless of slowly falling air temperatures, the upwarl 
and downward radiation currents become essentially 
constant. This further indicates a decrease in atmos 
pheric radiators. 


3. AVERAGE SOUNDINGS 
CLEAR WINTER 


An average cooling rate for clear winter ascents Wi 
obtained by averaging the upward, downward, and ne 
components of the radiative flux in 50-mb. steps for 2 
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Ficure 2.—Clear winter radiometer-sonde ascent. 


ascents. Cooling was then computed from these results 
providing a mean vertical distribution of the terrestrial 
radiation and its divergence. A simple computation, 
equation (1), gives the profile of cooling. In figure 4 this 
werage cooling is compared with Mdller’s [4] indirect 
estimate of the mean mid-latitude annual infrared cooling. 
Except for the observed decrease in cooling just beneath 
the mean height of the tropopause, the indirect estimates 
winecide favorably with the directly measured cooling. 
Iondon’s [5] computation of the mean March infrared 
‘ooling also lies well within a comparable range but is offset 
season. Between 150 and 200 mb. there is a tendency 
for constancy in RY and RA, though it is not as evident 
isin the average summer profiles. 


February 17, 1958, 2136 cst, Madison, Wis. 


CLEAR SUMMER 


The 12 summer flights, reaching a mean height of 75 mb., 
furnish a representative mid-latitude average of the verti- 
cal distribution of radiative fluxes and the infrared 
atmospheric cooling for a continental climate. The data 
of figure 5 are presented in the manner of figure 4. 

Between 70 and 150 mb. in the average summer sound- 
ings, the gradient of the downward- and upward-directed 
radiation fluxes tends to zero. This indicates that under 
clear summer conditions there exists a marked decrease 
in atmospheric radiators near the 100-mb. level and above. 
The mean cooling averages 1.7° C. per day compared with a 
mean winter cooling of 1.0° C. per day. 
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Figure 3.—Cloudy winter radiometer-sonde ascent. March 18, 1958, 2033 cst, Madison, Wis. Triangles represent Elsasser cooling 
evaluation. Circles indicate direct evaluation of radiative cooling. 


4. TERRESTRIAL RADIATION 


The values of this outgoing terrestrial (net) radiation 
for only the highest flights, clear and cloudy, in which the 
observed vertical gradient of the net radiation becomes 
zero, are averaged in table 1. For comparison, the esti- 
mates of Robinson [6] at Kew, Lettau [7] for Northern 
Hemisphere annual average, and Houghton [8] for North- 
ern Hemisphere annual average are included. 


TABLE 1.—Oulgoing terrestrial (net) radiation, langleys per minute 


Summer Annual Winter 
Madison 0.310 | Lettau------ 0.320 | Madison Average___ 0. 220 
Robinson (Kew) . 320 | Houghton. . 320 | Robinson (Kew)____ 270 
265 


5. CLOUD EFFECTS 


A characteristic phenomenon associated with ascents 
through clouds appears in figure 3 at 725 mb. Very strong 
radiative divergence immediately above a cloud top 
appears as an excellent indicator of the cloud top. This 
feature is obviously caused by the rapid decrease in the 
downward current of radiation measured by the upwaré- 
facing surface of the radiometer as it breaks into the 
clear night sky above the clouds. 


6. ELSASSER CHART COMPARISON 


The cooling profile indicated by triangular points in 
figure 3 represents an Elsasser evaluation for this flight 
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Figure 4.—Winter mean clear nocturnal ascent. 


Omitting that portion of the Elsasser [9] computation 
which is dependent upon assumptions made for the cloud 
emissivity, we find the general agreement between the 
direct and indirect cooling evaluation good. However, the 
humidity element was out of range above 450 mb., eliminat- 
ing an Elsasser evaluation for nearly two-thirds of the 
ascent. On the basis of the indirect chart radiation evalua- 
tions for a limited number of flights, the radiation values 
appear somewhat underestimated when compared with our 
direct ineasurements. This will be the subject of a forth- 
coming report. The magnitude of the directly measured 
warming and cooling below and above cloud bases and tops 
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Squares indicate annual cooling rate estimate by Moller. 


agrees with indirect estimates of Mller [4]. 


7. PRACTICAL APPLICATIONS 


Some of the apparent practical applications of the radi- 
ometer-sonde ascents are for cloud top determinations and 
in studies of cloud dissipation, propagation to the surface 
of high-level heating and cooling, stability changes, and 
the role of long-wave radiational destabilization as a 
trigger mechanism for nocturnal thunderstorms. The 


actual measurement of the infrared components under a 
variety of situations is providing a valuable tool with 
which these applications can be studied. 
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Figure 5.—Summer mean clear nocturnal ascent. 


8. SATELLITE METEOROLOGY IMPLICATIONS 


There is, of course, a natural tie between satellite radia- 
tion measurements and similar measurements with the 
radiometer-sonde. One might state that the balloon-borne 
radiometer fills in details integrated out by satellite 
measurements. Primarily, a radiometer-sonde network 
could augment the satellite measurements. Of necessity 
the satellite horizontally integrates radiation measure- 
ments of the earth and atmosphere over large areas because 
of its great speed and height. On the other hand, radi- 
ometer-sonde observations, though limited in height, pro- 
vide exceptional vertical detail over a limited area. The 
satellite measures the large-scale radiation balance of the 
earth-atmosphere system from a point in space. The radi- 
ometer-sonde measures radiative fluxes within the earth- 
atmosphere system. 


9. ACCURACY OF MEASUREMENTS 


Standard calibration procedure includes a separate base 
line check at three temperatures of the three thermistor 
used in each flight, together with a surface check of the 
actual net radiation against a suitable standard. 

Many additional tests and precautions, their actual de 
tails too extensive for this report, are either routinely 
performed in preparation for each flight or have been 
conducted as controlled experiments. Among these are 
the following: 


1. Dry nitrogen purging of each radiometer prior t_ 
flight to prevent internal moisture condensation. 

2. Individual calibration of the three thermistor sensor 
to prevent systematic errors. 

3. Use of proportional changes in both vertical and 
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horizontal dimensions of the radiometer as a check on 
yalidiiy of the conduction terms.. 

4. l'se of discrete changes in length of flight train as a 
check on solid angle intercept of balloon in affecting 
radiation measurement. 

5. Substitution of bead thermistors for standard rod 
thermistors to improve lag characteristics of radiation 
sensing surfaces. 

6. Extensive computations of transmissivity of poly- 
ethylene ventilation shield cells covering range from 0.8 
to 40 microns. 

7. Extending emissivity calculations for blackening 
agent, black “Decoret” Enamel, to 40 microns. 

8. Comparison, in flight, of single- and double-layer 
polyethylene ventilation shields as to their respective ef- 
fects on the ratio of conduction to convection transfers in 
these horizontal layers. 

9. Tests of effects of nonhorizontal orientation of radio- 
meter in flight by means of pendulum effect. 

10. Evaluation of validity of internal conduction term, 
(;,in equation (1) through its elimination in a comparison 
flight of a “disc” radiometer. 

11. Determination of the profile of the gradient of 
temperature in the horizontal through the radiometer 
through multiple sensors. 

The effects of a flight trajectory through clouds on the 
absolute accuracy of either 2Y or LA have not been com- 
pletely determined. While it is believed that such effects 
on radiation measurements as moisture on the top poly- 
ethylene shield are rapidly eliminated due to ventilation, 
this is being carefully investigated. The rapid decrease 
inthe magnitude of the downward propagating radiation 
current indicates a diminution in this assumed adverse 
effect immediately upon exit from cloud. 

A careful error analysis has been made of the accuracy 
of temperature measurement including these random 
errors: thermistor error, thermistor lag, recorder error, 
calibration error, and transmitter error [10]. The “ac- 
curacy” ascribed to any measuring instrument may be 
misleading because of the many ways of estimating and 
expressing accuracy. Customarily, we employ various 
statistical measures of dispersion to indicate this accuracy 
ind those most often used are the standard deviation, the 
iverage deviation, and the probable error. Their inter- 
relation is well known. However, we will limit our error 
discussion to the standard deviation, o. The root of the 
sum of squares of individually computed temperature 
measurement standard deviations in the error analysis 
mentioned above provides a value for one standard devia- 
tion, o, of 0.4° C. Setting out confidence limits at 95 
percent, we must accept a maximum error (20) of 0.8° C. 
Such x random error becomes most important in comput- 
ing atinospheric cooling and can possibly result in maxi- 
um cooling errors of up to +1.2° C. This maximum 
cooliny rate error is large in comparison to some of the 
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observed mean values, but it should be remembered that 
95 percent of the errors will be less than 1.2° C. and 67 
percent less than 0.6° C. We are experimenting with 
methods of reducing this maximum error to +0.5° C. 
Such methods center primarily about a measurement of 
AT’, the temperature difference between top and bottom 
sensors, rather than the absolute top and bottom tempera- 
tures, which can go far toward eliminating the inherent 


error in the small difference between two large numbers. 


Consider the expression for net radiation, ?y, 
Ry=o(T}—T}) 


+1.37 | (2) 


Investigation reveals the relative size of the first and 
second terms on the right, the bleck-body and conduction 
terms, to be in the ratio 1.0 to 0.6 in the lower flight levels. 
At the tropopause this ratio changes to 1.0 to 1.8. Above 
the tropopause with the vertical change in the upward 
and downward radiation currents tending to zero, the lag 
term, third on right, becomes insignificant. The maximum 
size of the lag term never exceeds 5 percent of the net 
radiation and this occurs only in cases of extreme gradient 
of the observed top and bottom radiating temperatures. 
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THE DISTRIBUTIONS OF FREEZE-DATE AND FREEZE-FREE PERIOD 
FOR CLIMATOLOGICAL SERIES WITH FREEZELESS YEARS* 
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(Manuscript received December 3, 1958; revised March 4, 1959] 


ABSTRACT 


It is shown that the freeze distribution is a mixture of the distribution of freeze-date and the simple dichoto- 


mous distribution of freeze and freezeless years. 


This is applied both nonparametrically and assuming a normal 


distribution of freeze date to three stations at three different thresholds to obtain the probabilities of freeze 
before or after any date. The distribution of the freeze-free period is developed and application made to one of 


the stations to obtain probabilities of the freeze-free period being less than a given time interval. 


The expressions 


for the mean freeze-date and freeze-free period are also developed and estimates made for the stations treated. 


1. INTRODUCTION 


The estimation of freeze probabilities from complete 
freeze-date series has been treated by Thom and Shaw 
[1]. When the freeze-date series for an observation sta- 
tion is éncomplete in the sense that some years experienced 
no freeze, there is, of course, a probability of freezeless 
years. This is common in more southerly latitudes, espe- 
cially for freeze thresholds below 32° F. With the addi- 
tion of the no-freeze probability component, a quite dif- 
ferent problem in the estimation of freeze probability 
arises. This has been discussed by Spillman et al. [2]. 
They gave rules for finding the mean recurrence interval 
for the incomplete series; but since they did not recognize 
the more general statistical aspects of the problem, their 
rules are not completely convertible to probability state- 
ments. 


2. THE FREEZE DISTRIBUTION FUNCTION 


The model for determining freeze probability may be 
thought of as a mixture of two distributions: one a dis- 
crete distribution of no-freeze and freeze, the other an 
essentially continuous distribution of freeze-date for years 
when freeze occurred. In this discussion, the period over 
which spring freeze-date is assumed to range is January 
1 to June 30, and that for fall freeze is from July 1 to 
December 31. These are arbitrary, and other dates may 
be assumed if it suits a particular purpose better, as we 
shall see later. The model is seen to be equivalent to con- 
centrating a probability of no-fgeeze at an arbitrary point 
before the beginning of the season for spring freeze and 
after the season for fall freeze. 

We define the spring freeze-date series as in [1] to be 
the series of annual last dates in spring on which a min- 
imum temperature less than the threshold temperature 


*This paper is based on work done while the writer was Visiting Pro- 
fessor of Statistics, Biometrics Unit, Cornell University [3]. 


(32°, 28°, 24°, 20°, 16°) has occurred. The fall freeze- 
date series is defined by substituting the words fall for 
spring and first date for last date. . 

On the basis of the results given in [1], we shall assume 
that the climatological series comprising freeze-dates 
mixed with no-freeze occurrences are random variables. 
It follows then that the distribution functions of freeze 
may be found, and that these will completely define the 
freeze series populations. 

The distribution function is defined as usual by 


F(2)= fu) du (1) 


where f(w) is the probability density function (pdf) or 
frequency distribution and F(—o)=0 and F(«)=1 
Here/ (x) is the probability that wu is less than @; and 
when « is continuous this is identical with the probability 
of a value less than or equal to # Clearly in the spring 
we shall be most interested in the probability of a freeze 
occurring after a, and hence, we shall be interested in the 
form 1—F (a) which gives that probability. In fall we 
shall be interested in the probability of a freeze before # 
which is given by F (a) itself. 

The mixed distribution of freeze-date and no-freeze for 
spring may now be derived as follows: Let g, be the prob- 
ability of no-freeze occurring in spring according to the 
model assumed above. 1—g.=pz is then the probability 
of a freeze after the beginning of the freeze season. Ac 
cording to the definition of mixed distributions [4], the 
distribution function for spring freeze will be 


= gst ps F's(a) (2) 


where /’s(a:) is the distribution function of spring-freez 
date when freeze occurred. We employ large S and A te 
indicate the continuous portion of #, and small a and $ 
to indicate the discrete portion. It is seen that G(7) 8 
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a distribution function, for if # takes a small value 6 by 
definition /(b)~0, meaning that no freeze-date will 
occur before 6 where G(b)=q.. This is the probability 
that no-freeze will occur. If a is not a member of either 
sor 8, then G(a)=0, for it is impossible that either one 
of the events freeze or no-freeze should occur. On the 
other hand, if # is a large value, c, then by definition 
F(c) ~1, and which is the probability 
that either no-freeze or a freeze-date have occurred before 
¢, Thus (2) is a distribution function. It will be clear 
now that the probability of a freeze after date @ will 
depend not only on 1—F'(x), the probability of freeze- 
date when freeze has occurred, but also on the probability 
that a freeze will occur at all on any date. This, of 
course, is 1—gs Or Ps. 

As we have seen, G(a) of equation (2) gives the prob- 
ability of freeze or no-freeze before x, whereas our main 
interest is in the probability of freezes after x. This is 
clearly one minus the probability obtained from (2). Let 


H(x)=1—G(z), 
then 


H(z) —qs— pF, (2) ; 


and since p,+q,=1, 
H(z) =p,|1— Fs(z)). (3) 


If we write 


I(x) =1—F;(z) (4) 
equation (3) becomes 
H (x) (2). (5) 


This gives the probability of a freeze occurring, after 
date x in spring. 

For fall freeze we have a similar mixed distribution 
except that the probability of no-freeze is now concen- 
trated after the fall freeze season and again does not 
enter into the probability before #. Hence the date 
distribution is 


J (x) = pala (2) (6) 


Here p, is the probability of a fall freeze and /'4(#) is 
the distribution function on date. Since the distribution 
function gives the probability of freeze before date 2, 
equation (6) gives the required probability directly. If 
the probability after v is needed, this may be obtained 
from 1—J(a#). This then includes q., the probability of 
no-freeze in fall (autumn). 

It should be noted that equations (5) and (6) hold 
generally, for in the situation where freeze occurs every 
year, as discussed in [1], g=0 and p=1. 


3. ESTIMATION OF FREEZE PROBABILITIES 


The main objective in developing the freeze distribu- 
tion is to provide the means of obtaining probabilities. 
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Thus, proper estimation of the terms in (5) and (6) will 
provide estimates of the required probabilities. There 
are two ways in which we can estimate these terms: Hav- 
ing estimated the p’s, we may estimate the / and /’, fune- 
tions directly from the data, or we may first estimate the 
parameters of the 7 and F, functions. Before we can 
perform these estimations, however, we must define our 
climatological variable, freeze-date, more closely. 

Clearly calendar date would be unsatisfactory as a 
variable. However, we can easily convert calendar date 
to day number beginning from some suitable base date. 
This will facilitate computations, the statistics from 
which may be readily converted back to calendar date. 
Since freeze dates vary over the periods July 1 to De- 
cember 31 and January 1 to June 30, January 1 has been 
chosen as the base date. In leap years the 366-day year 
was employed. 

Inasmuch as the base date will affect the mean of an in- 
complete freeze series, it might appear to be somewhat 
better to place the base date at a point halfway between 
the means of the fall and spring dates when freeze actually 
occurred. However, this would result in little refinement 
and would cause great inconvenience, for a computation of 
the halfway date would be required for each station. 
An examination of a number of stations showed that the 
halfway date usually occurs a few days after January 1. 
In view of the larger dispersion of the spring dates, the 
ideal base date, on probability considerations, should be 
displaced backward in time somewhat from the halfway 
date. This, together with the fact that the choice of base 
date does not greatly affect the probabilities, seemed to 
make the January 1 base the most satisfactory. All data 
with which we shall be concerned have therefore been 
coded to January 1. Tables 2, 3, and 4 show the freeze 
dates coded in this manner for Anniston, Birmingham, 
and Auburn, Ala. The data are arranged in order of 
increasing date and labeled with order number %. The 
three southern stations were chosen to emphasize the in- 
completeness aspect of the freeze series which is the central 
problem of the present analysis. } 

Our first estimates of 7 and J will be empirical or non- 
parametric. These involve first the estimation of / and F',. 
While ordinarily with quite long series these would be 
estimated by 4/m where & is the order number and m is 
the number of actual freeze dates, it has been found that 
for a continuous distribution the following equation gives 
estimates which are more unbiased at the smaller and 
larger probabilities : 


(7) 
For spring this becomes 
k 
(8) 
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TaBLe 1.—Statistics for various freeze thresholds for Anniston, 
Auburn, and Birmingham, Ala. 
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TaBLe 2.—Estimated probabilities of 16° freeze for Anniston afte 
given dates in spring and before given dates in fall 


Freeze n | m and | 


(code) (date) | (date) 
ANNISTON 
Spring 
32 2) 8/80) 14.0 | 1.000 | 0.820 |  —0.195 
23) 2» 29/ 70.1) 3/1!) 3/11) 169) 1.000! .768 005 
24 29; 28) 55.0) 2/24| 2722) 180| .966| .779) —.846 
20 29 47.0) 2/1| 192] | | 375 
16 16) 39.8) 29) 122) 18.8) .552| .736 431 
Fall 
32, 28 28 310.3) 11/6) 11/6 | 1.000 | 0. 795 0.310 
2s 28 28 323.8) 11/20) 11/20) 122) 1.000)  .770 —. 985 
24 2s 22) 333.4) 11/29| 126) 139| .786| 888 
20 2s 14) 12/78!) 10.7) .500 | 829 | .279 
16 2s 9) 345.1 12/11 | 12/25) 117) —.710 
AUBURN 
Spring 
32; 30 30 | 80.3 3/21) 3/21 | 16.6 | 1.000 | 0.805| —o.219 
28 30; 181] 1.000| —.478 
24 2/14) 24) .981] .813 | —. 381 
20 2; 426! #212; 23| 196| .864 —.173 
12] 320 wis} 151) | 789 | 
Fall 
32; 317.9; 11/6 | 9.6 | 0. 967 | 0.42 | 
28 30 27 | 333.1 | 11/20) 12/2| 11.4] .900| .725| 910 
28 30 19 | 338.6) 12/5| 12/4) | 468 
20 30 | 10; 340.1 | 12/6) 12/23 1.3) 807 | —. 194 
16 30 6| 12/13 | 1227 | 113) 738 —. 266 
| | 
BIRMINGHAM 
Spring 
$2; 77.6 | 3/19 | 3/19 | 16.3 | 1.000| 0.764 0.493 
2s 30} 627| 3/4| 34| 17.3] 1.000] .770| .097 
24 30) 50.5 | 2/20) | 15.7] | BIS 
20 30; 19| 1/28| 17.0] :633| [788 —. 219 
16 30 4 | 34.4 | 2/3 | 1/16 | 19.7| .467| .849 170 
Fall 
32 | 30 29 | 316.7 | | 12.2) 0.967 | 0.797) —0.019 
23| 30 28 | 334.7/ 12/1 14.3) 825 116 
30 21 | 339.8/ 12/6) 1213| 125] .700| .854 362 
30 12} 3421 12/8) 1272) (400! —. 531 
16) 8 | 346.5 | 12/13 | 12/26 | 11.5) .267| .819 —. 135 


The star indicates a nonparametric or distribution-free 
estimate of a parameter from a sample. 

To complete the estimation of 7 and J, we must esti- 
mate p, and p. Since the freeze, no-freeze series forms 
a discrete distribution the estimates are found from 


(9) 
where p is the parametric estimate of p,m is the number of 
years with freeze, and 7 is the number of years with freeze 
or no-freeze. The number of years with no-freeze is, of 
course, 

The estimates p, and p, are shown in the p column of 
table 1. These were obtained by applying equation (9) 
to the m’s and n’s listed there. The statistics for all thresh- 
olds are given in table 1, although only Anniston 16°, 
Auburn 24°, and Birmingham 20° are discussed in full. 


Spring 
k | Date r ye | I H* fr 

1/7 | 7 0.941 | 0. 959 0. 519 0. 529 
1/16 | 16 824 | 896 495 
1/25 | 25 765 | 782 | . 422 
1/28 | 28 | | 732 | 390 ‘404 
1/31 31 647 | 357 ‘3% 
8. 2/10 | 41) 472 | 
9.. 2/11 | 42 | 471 452 260 "250 
| 44 412 409 | 227 2h 
2/19 | 50 . 294 | 291 . 162 
2/27 | 58 235 | 166 . 130 “092 
14 2/28 59 176 | . 154 097 085 
15_. 3/4 63 118 | 107 | 065 058 
3/10 | 70 | 059 | 054 033 030 

Max !d|=0.116 
Fall 
k | Date r | J 

1. 11/16 321) 0.100} 0.018 | 0.082 0.006 
12/ 3 337 | 200 | 251 064 
12/ 7 341 300 | 367 026 118 
4. 12/9 343 400 433 | 128 138 
5.. 316 500 | 532 | . 161 
351 700 | 691 225 2 
12/23 | 357 | 800 | 844 257 271 
12/26 B | 


Max |d|=0.082 


The 7* and /'4* are estimated by applying equations (7) 
and (8) to the %’s and m’s of tables 2, 3, and 4, giving the /* 
and #* columns of spring and fall freeze of those tables. 
From equations (5) and (6) it is seen that it is necessary 
to multiply the 7/* and F* by p, and pu, respectively, to 
obtain H* and /*, the nonparametric estimates of the 
mixed distribution. From table 1 we find for Anniston 
16°, ps=0.552 and pa=0.321. Multiplying these, respec- 
tively, by the values of 7* and F* from table 2 gives the 
H* and J* columns of the table. A similar calculation 
applies to tables 3 and 4. It is seen then that /7* gives the 
probability that a freeze occurs after date # in spring 
and /* the probability that a freeze occurs before date # 
in the fall. From table 2 we see that the probability ofa 
16° freeze at Anniston after February 11 is 0.260, after 
March 10 it is only 0.033, or about 1 year in 30. The 
probability of a 24° freeze occurring before November 19 
at Auburn is 0.063 from table 3, and of a 20° freeze after 
March 10 at Birmingham is also 0.063 from table 4. We 
do not recommend the nonparametric estimates for use in 
obtaining probabilities since we have a theoretical distribu- 
tion as we shall see below. However, the empirical proba- 
bilities are necessary for judging the fit of the theoretical 
distribution so they are plotted in distribution function 
form as broken lines in figures 1, 2, and 3. If one desires 
to use the empirical or nonparametric probabilities, the 


_recommended form of graph would be that shown in these 


figures. 


‘ 
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TasLe ‘Estimated probabilities of 24° freeze at Auburn after given 
dates in spring or before given da'es in fall 
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TABLE 4.—Estimated probabilities of 20° freeze at Birmingham after 
given dates in spring or before given dates in fall 


Spring Spring 
k Date i H* k Date r Ir H* 
1 1 0. 964 0. 980 0. 897 0. 912 a ee BS 1/15 15 0. 950 0. 970 0. 601 0.614 
1/ 2 2 . 929 .977 . 865 . 910 EER eee 1/16 16 . 900 . 952 . 570 . 608 
1/15 15 . 893 918 831 855 tons 1/24 24 . 850 . 875 . 538 
1/17 17 . 857 . 903 . 798 . 841 1/28 28 800 821 . 506 . 520 
1/24 24 821 834 . 764 . 776 1/30 30 . 750 . 791 . 
1/28 28 . 7386 . 785 . 732 . 731 A See 1/31 31 . 700 .773 . 443 . 489 
2/1 32 . 750 . 726 . 698 . 676 2/10 41 . 650 . 571 411 61 
2/1 32 .714 . 726 . 665 . 676 2/1 42 . 600 . 48 B80 47 
2/ 3 34 . 679 . 695 . 6832 647 2/13 4d . 550 . . . 319 
2/8 39 . 643 . 06 . 599 . 564 2/15 46 . 500 . 456 .317 239 
2/9 40 007 . 1 . 565 . 550 2/19 50 . 450 . 367 285 . 232 
2/11 42 . 571 . 552 . 532 514 ae 2/19 50 . 400 . 367 . 253 . 232 
2/12 43 . 536 . 536 . 499 499 Sa a 2/19 5D . 350 . 67 . 222 . 232 
2/13 44 . 500 . 516 . 466 . 480 2/19 . 300 067 . 190 . 232 
2/19 50 . 464 . 405 . 432 .377 St ee 2/28 59 . 250 . 198 . 158 . 125 
2/19 50 429 . 405 . 399 .377 3/3 62 . 200 . 14 . 097 
2/23 54 . 393 334 . 366 3/4 63 150 .140 005 089 
2/28 59 . 357 . 255 . 332 . 237 Be 3/10 70 .100 . 076 . 063 048 
2/28 59 . 321 . 255 . 299 . 237 Dittdducdcndadaducdieba 3/14 73 050 049 . 032 . 031 
3/1 60 . 236 . 239 . 266 . 223 
3/8 67 "179 "149 "167 "139 Max |d|=0.083 
3/11 70 . 148 . 133 lll 
3/13 73 .107 .102 . 100 095 Fall 
3/14 73 . O71 . 093 . 066 . 087 
3/20 79 036 | 055 034 051 
k Date z Fr F J* J 
Max/d|—0.102 
Re eee 11/15 320 0.077 0. 021 0. 031 0. 0O8 
11/24 328 . 14 . 109 . 062 
Fall 11/29 333 | . 215 . 092 . 086 
12/ 2 337 . 308 . 802 .121 
12/ 6 340 . 433 . 14 .173 
k Date Ft J* J 12/ 8 342 462 "504 185 ‘202 
12/11 345 . . 606 . 215 . 242 
12/15 349 615 . 736 246 204 
11/15 320 0. 050 0.077 0. 032 0.049 12/16 350 . 692 . 764 .277 . 306 
11/19 323 . 100 . 063 074 12/16 350 . 769 764 308 . 306 
11/20 324 . 150 . 133 095 OS4 12/21 355 46 877 338 . B51 
11/24 328 . 200 . 212 134 12/22 356 . 923 sud 369 358 
11/25 329 . 250 . 233 . 158 147 
11/26 331 . 300 . 258 . 199 163 
11/27 331 . 350 . 234 . 222 180 Max =0.121 
11/28 332 . 400 . 253 197 
11/29 333 . 450 . 337 . 285 213 
11/30 334 . 500 . 367 .317 232 
eae cee 12/ 2 337 . 550 . 429 . 348 272 
12/ 9 343 . 600 . 641 380 406 
40 continuous component of the mixed distribution for nor- 
= 517 —mality which is fitted to the date of freeze in the series of 
1400 ‘oor | , ces actual freeze occurrences. 
Probability tables used in [1] due to Geary [6] are 


Reed [5] and later Thom and Shaw [1] found that the 
normal distribution provided very good fits to freeze- 
date series under a wide range of conditions for the 32° 
and other thresholds for complete series; i.e., for p=1. 
The series, of course, tend to be more incomplete the far- 
ther south we go. This is also accompanied by a shift in 
the center of the distribution toward the colder season; 
ie, toward winter from both fall and spring. This shift 
of the distribution center naturally causes some concern 
since the tails of the distributions on the winter end could 
begin to show the effect of boundedness, and hence depar- 
ture from normality due to the decrease in time interval 
over which late fall and early spring freeze can range. 
This was also the reason for testing our theory on stations 
ina southern region where conditions are most stringent. 
To verify a part of our theory it is necessary to test the 


»08099—59——4 


again employed to test for normality. In these, a, the 
standardized mean absolute deviation from the mean, and 
\/b,, the standardized central third moment, are measures 
of kurtosis and skewness. These statistics are listed in 
table 1. Using the tables of [6] it was found that none 
of the a’s are significant at the 0.10 probability level, and 
only the \/d, (in italics in table 1) for 28° freeze in fall 
at Anniston is significant at the 0.02 level. The four 
largest values of \/d, are individually significant at the 
0.10 level but average near zero. Two of these are nega- 
tive and two are positive; however, fall and spring each 
have a negative and positive value. This is in disagree- 
ment with what we would expect on the basis of the pos- 
sible boundedness mentioned above which would cause 
negative skewness in fall and positive skewness in spring. 
We feel, therefore, that it is reasonable to assume that 
these larger values were a result of sampling and that, 
therefore, the normal distribution satisfactorily fits the 
continuous component of the mixed distribution of freeze 
occurrence and date. 
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10 7 gives the probability of freeze occurring after date g jy fi 
ANNISTON spring and the probability of an occurrence befop 
ot plicit rs. | date a in fall, both on the condition that freeze has actually [Bg 
occurred. These probabilities are parametric estimates qe 
rs — Paremetvte : and are indicated by / and F in tables 2, 3, and 4. The di 
@ normal estimates and the nonparametric estimates of the al 
ist | probabilities may be compared by contrasting the careted p 
and the starred distributions. 
: SPRING 4 Although there is little question about the adequacy of in 
4 the normal distribution in fitting freeze-date, it is of in. pe 
terest to test the fit in another manner. For this purpose in 
8/80 we use the easily applied Kolmogorov-Smirnov test for bn 
DATE which Massey has provided convenient tables. Massey 3 
Ficure 1.—Probabilities of 16° freeze at Anniston, Ala., occurring [7] has also examined the power of the test and has found 
’ after any given date in spring or before any given date in fall. it superior to the x? test in the cases analyzed. The test be 
is carried out by examining the significance of dk 
sei 
T T T max | d| =max | k/m—ZJ(@) |. (10) an 
BIRMINGHAM sa 
ak 20° FREEZE : This is the maximum absolute difference between /:/m and in 
the normal distribution function. The maxima of |d| no 
are shown at the foot of each distribution table. Each sp 
ng 7 value may, of course, be tested separately; however, we 
Z may test them all at once by considering the maximum of ou 
© 4 | the max |d| for the longest record employed in fitting tot 
the normal distributions. This will be a more stringent Ls 
; ? | _ test than if we had used the actual length of record, length W 
7 m. Max (max |d|) we see to be 0.133 and the longest an 
record is 19 years. For these arguments Massey’s table dis 
1730 2/19 vio gives (1 
DATE P(max | d| >0.133) >0.20 fol 
Ficure 2.—Probabilities of 20° freeze at Birmingham, Ala., occur- an 
ring after any given date in spring or before any given date in fall. Since 0.20 is a rather large probability, the fit is good; in (0 
fact, all the fits are as good or better than this. This the 
strengthens our conclusion of normality reached above. 
10 . The goodness of fit also extends to the mixed distributions 
AUBURN since there is little question of the fit of the p’s. 
Pm eo FREEZE : The mixed distributions of spring and fall freeze are = 
obtained from equations (5) and (6). These are the 1's 
ee and J's of tables 2, 3, and 4 and are obtained from the Is 
and F’s by multiplying respectively, by ps and Th 
8 | pa | and J* are the nonparametric mixed distributions, and int 
e; H and / are the parametric mixed distributions. These itie 
are plotted in figures 1, 2, and 3. Here one may observe dis 
we the rather good fits of the smooth mixed theoretical dis- (36 
4 tribution to the broken line empirical distributions. 
Probabilities of freeze before or after any date may be 
wit, read from the smooth curves. 
Figure 3.—Probabilities of 24° freeze at Auburn, Ala., occurring and 
after any given date in spring or before any given date in fall. 4. MEAN FREEZE DATE 
Although the mean freeze date of the incomplete freeze 
I(v) and F(«#) were fitted as normal distributions to series is more difficult to interpret, it is perhaps of some Her 
the spring and fall freeze-date series in the usual manner formal importance to consider it. Ordinarily the mean don 


by estimating the means and standard deviations. Z(#) value of even a mixed distribution is obtained readily by don 
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finding the expected value of the distribution. Here, how- 
ever, (he expected value depends on what value we assign 
to the dates on which we assume the probabilities g, and 
q, to be concentrated. We made some provision for this 
difficulty when we considered the freeze distributions 
above. but this was not as important there because the 
probabilities were not affected by the base date. If we 
make the reasonable assumption that the maximum grow- 
ing season is 365 days, then it follows that when no-freeze 
occurs in spring, but one occurs at date a in fall, the grow- 
ing season is # days long. Ifa freeze occurs at x in spring, 
but none occurs in fall, then the growing season is (365 
-z) days long. We have previously chosen our day num- 
ber code to give a reasonable interpretation to the mean 
dates; hence if we add our assumptions about growing 
gason length, g. will be assumed concentrated at code 0 
and g, at code (365+0). 0 and (365+0) are actually the 
same hypothetical day chosen so as to fit the conditions 
imposed by our model of the freeze-free season. We may 
now readily define the mean or expected values of the 
spring and fall freeze date. 

For convenience we designate as and a, as the continu- 
ous parts of the spring and fall freeze variables on the 
total interval (1,365) and a and a, as the discrete parts. 
z, only takes the value (365+0) and a, only the value 0. 
We define the pdf's on data in spring and fall to be 7 (as) 
and g(w4), respectively. These are the derivatives of 
distribution functions expressed generally by equation 
(1). We need also to define the mixed variable date u 
for spring and v for fall, so as to include both the discrete 
and continuous parts. wu and v then have the total range 
(0, 365+0). With these additions, we may now express 
the mixed pdf of spring freeze by 


h(u) = st ps f (ws) (11) 
and for fall freeze by 
j(v) = Pag + qu (12) 


The expected value or mean is defined as the sum and 
integral of the products of the variates by their probabil- 
ities. Taking into account the fact that g, and ga are 
discrete components of probability concentrated at 0 and 
(365+0), we find the mean value for spring to be 


E(u)=0Xq.+p; rsf (xs)drg (13) 
and for fall 


E(v) = (14) 


Here \ is the domain of spring freeze pdf and A is the 
domain of fall freeze pdf. Although in reality the 
domain (S+A) is (1, 365) as mentioned previously, we 
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find it convenient to assume S and A to be of infinite 
extent. This is plausible because the probability densities 
become very small near #=1, 365, and 364/2. With this 
assumption, together with the previous one that f and g 
are normal pdf's, it is clear that the integrals in (13) and 
(14) are mean values conditional on the occurrence of 
freeze. These may be expressed by (as) and F(a,) 
meaning the expected values in domains S and A. Equa- 
tions (13) and (14) then reduce to 


E(u) = ps E (as) (15) 


and 


E =365 gat poll (a1). (16) 
Inasmuch as /'(#s) and /'(a,4) are the means of normal 
distributions, they will be best estimated by the arith- 
metic means of spring and fall freeze dates ay and a4. 
Substituting these estimates together with the estimates 
of the p’s and q’s, we find the estimated mixed means for 


spring 


u= sy (17) 


and for fall 


V=365 gat Pals. (18) 
These are listed in table 1 for the stations studied and 
have been tabulated for a large number of Weather 
Bureau stations. 


5. FREEZE-FREE SEASON DISTRIBUTION 


It will be clear from our definition of the freeze-date 
variable for spring and fall that we may express the 
freeze-free season by the variable y=v—wu. We shall 
assume, based on previous work for complete series {1}, 
that spring and fall freeze dates are independently dis- 
tributed. This seems all the more justified for incom- 
plete freeze series since the spring and fall distributions 
are located farther apart in time than those for complete 
series. In order to find the distribution of freeze-free 
season, we must obtain the distribution of y. This may 
be most conveniently done by moment-generating func- 
tions or characteristic functions (cf). The former are 
ordinarily Laplace transforms and the latter Fourier 
transforms. It makes little difference which we use here 
since our development need only be a synthesis of known 
results. We prefer the cf because it has somewhat more 
general application. 

The cf of a pdf is defined as the expected value of the 
Fourier kernel. (See [8], ch. 10.) 


o(t)= Ele) = et (19) 
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Here ¢ is an arbitrary variable and f(z) is the pdf. 
The analytical power of the ef arises from the fact that 
g(t) is a Fourier transform, and therefore the inverse 
may immediately be expressed as 


which gives the required pdf. The cf’s are obtained 
from the pdf’s of equation (11) and (12) (see [8], pp. 
57-58) by multiplying by the Fourier kernel and taking 


expected values. This gives 

p(t) = dst ps vs(t) (21) 
and 

w(t) =gat po pa(t) (22) 


For convenience in convolution we may take y=v 
+(-—w) in which case the variable —w has the cf 
ys(—t) ([8], p. 185), and we may rewrite (21) as 


o(t)=qstps ps(—t). 
By the principle of convolution of distributions the cf of 
the sum y=v+(—w) will be the product of (22) and 
(23) ;1.e., 


Q(t) =Geqat + qapeps(—t) 

+ pspa pa(t)ps(—t). (24) 
Assuming normal distributions for spring and fall freeze- 
date, we consider the last term first. This is the term 
involving freeze occurrence in both spring and fall. It is 
well known that the cf of a normal distribution with 
mean » and standard deviation o is 


w(t) =exp | (25) 
using our variable for fall freeze. Hence for spring we 
have the required cf 


é(t) =exp og” — | (26) 
Multiplying (25) by (26), we find the cf of the mixture 
component of freeze-date to be 
of=exp [ — (04? +05") +it(ps—ps) (27) 
Since this is of the form (25) with mean ps—ps and 
variance (o47+ 3"), it is a normal distribution with these 
parameters. We shall express the pdf of this cf as w(z) 
where z= (@4—@s). 

As to the other terms on the right of equation (24), 
the first term is the product of the probabilities that there 
is no freeze in either spring or fall and is consequently 
the probability that the freeze-free season is 365 days. The 
pdf's of the second and third terms may be interpreted 
as follows: The second term is the mixture component 
with no-freeze in spring but one in fall; hence the variable 
2=(a#4—0)=a,. Since the cf is for the pdf g, this remains 


(23) 
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g(wa). The third term of (24) is the mixture component 
with freeze in spring but none in fall. For this the freeze. 
free season variable is 2=365—as and the pdf becomes 
f(365—ag) which has the cf gs(—¢) with a location shift 
to account for our chosen freeze-date scale. 

With these interpretations made, we may now write the 
pdf of freeze-free period as 


UY) = (4) + qapef (365—as) + pspaw(2). (28) 
If we integrate equation (28) term by term, remembering 


that the first term is a constant, we find the distribution 
function of freeze-free period to be 


Q(y) + as) + papaW (2) + 
(29) 

As we have seen previously, the distribution function 
gives the probability of an occurrence of freeze before date 
w; hence Q(y) gives the probability of a growing season 
less than y. This may be obtained from equation (29) by 
substituting the estimates for the p’s, g’s, means (a), and 
standard deviations (s) from table 1. 

As an example, we obtain @(v) for the 16° —16° freeze- 
free period at Anniston. This cannot be expressed as a 
single function; it must be compiled by adding the com- 
ponent probabilities of equation (29). From table 1 we 
find the mean of the 16° fall-freeze variate to be %4=345.1 
and the standard deviation s4=11.7. According to our 
definition, the mean for the spring freeze variate is 
365—Z 5. From table 1, 16° spring freeze, we find 
#s=39.8; hence, 365—%#s;=325.2. Since the standard de 
viation is not affected by either the algebraic sign of ¢; 
or its subtraction from 365, ss is also found from table | 
and is 18.8. From equation (27) we see that the sample 
mean of the last component is 2=%4—#%s=305.3, and the 
standard deviation is 


(11.7)? + (18.8)?=22.1. 


From table 1 we also find p.=0.552 and pa=0.321. The 
q's are one minus the p’s; hence g,=0.448 and ga=0.67%. 
These readily yield the coefficients in (29): ¢spa=0.14, 
JaPs=0.375, PsPa=0.177. The probability of a 16°-16 
freeze-free season of 365 days is gsga=0.304. Note that the 
coefficients of the probability function (28) add to unity 
as they should. 

The probabilities for each component of the mixture for 
a convenient set of freeze-free season durations were com 
puted using the above data and are listed in table 5. The 
theoretical distribution Q is listed in the eighth colum 
and is the sum along rows of the three components. The 
empirical distribution Q@* computed from the original 
freeze-free season series by equations (7) and (29) #8 
shown in the last column for comparison purposes. 

The theoretical and empirical freeze-free distributions 
are shown in figure 4, the former by the smoot} curve 
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Tarte 5.—Distributions Q and Q*) of freeze-free season durations for 16° threshold for Anniston, Ala., and probabilities for each component 
of the mixture 


aA A a A A A A | | 
y W Fa Fs PePaW | | | Q | F* | 
| 
250 0. 006 0.001 | 0,001 287 0.048 | 0.034 
255 002 288 095 O68 
260 . 020 . 004 . 004 290 . 148 102 
265 . 034 0. 001 . 006 . 006 295 . 190 136 
270 055 0.001 301 170 
275 085 ‘O15 315 204 
280 025 315 333 238 
295 ‘054 056 324 476 340 
300 . 405 . 090 . 072 . 034 . 106 324 . 524 374 
305 . 496 . 142 . OBS . 053 .141 327 . 571 40s 
310 . 583 0.001 . 209 . 108 . 078 . 181 337 . 619 442 
315 . 670 . 005 . 295 119 0. 001 ill . 231 337 . 67 476 
325 043 144 186 336 341 76: 544 
330 . 869 . 099 . 603 . 154 . 226 304 348 . 810 579 
345 964 496 320 562 357 952 680 
350 978 . 663 . 907 .173 . 095 . 340 . 608 
355 UNS . 802 . 944 .175 115 . 354 . 644 | 
360 993 . 898 . 968 .176 .129 . 363 . 668 
365 997 983 138 369 683 


| 
| 


and the latter by the broken line. The 365-day component 
of probability is the vertical line at right end of the 
theoretical curve. The maximum absolute difference be- 
tween the empirical and theoretical curves, adjusted to 
unity to make the continuous part a distribution function, 
is 0.089. For sample size 20, the number years with 
freeze, Massey's [7] table gives a much larger value 0.231 
at the 0.20 probability limit. The fit of the theoretical 
distribution to the actual data is therefore very good. 
Probabilities that the 16°-16° freeze-free season is less 
than any number of days read on the abscissa may be 
read from the ordinate of the figure. 
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ficure 4.—F'reeze-free distributions (16° threshold) for Anniston, 
Ala. 


6. MEAN FREEZE-FREE SEASON 


The mean or expected value of the freeze-free season 
which is at least of some formal interest can be readily 
found by the usual methods from equation (28). Multi- 
plying each term by its variate and integrating, recalling 
that the first term is associated only with 365, and sub- 
stituting rough values, we find 


F=365 Fat Gaps (365—Bs) + (30) 


It was pointed out to me by Dan Harton that y can 
obviously be obtained much more easily from the difference 
of equations (15) and (16). This gives 


Y= 365 gat Pa La— Ps (31) 


By a considerable amount of algebraic manipulation, 
equation (30) may be reduced to equation (31). This is 
of interest since it shows that convolution has produced 
a distribution (28) consistent with the basic assumptions 
as, of course, it must if it is correctly defined. 

Using either formula together with the values found 
above, we find 7, the mean 16°-16° freeze-free season at 
Anniston, to be 336.6 days. The means of spring and fall 
freeze and freeze-free period are available from Weather 
Bureau State climatologists for a large number of stations 
in the United States. 
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THE LOWEST TEMPERATURE IN ANTARCTICA 


NINA A. STEPANOVA 
Office of Climatology, U.S. Weather Bureau, Washington, D.C. 
{Manuscript received April 17, 1959] 


The IGY period 1957-58 was the first time in history 
that continuous meteorological observations were made 
in the interior of Antarctica. The U.S. IGY Amundsen- 
Scott Station (90° S., elev. 2,800 m.) and the Soviet IGY 
stations, Vostok (78°27’ S., 106°52’ E., elev. 3,420 m.) 
and Sovetskaya (78°24’ S., 87°35’ E., elev. 3,570 m.) are 
located on a high plateau of the Antarctic Continent and 
the minimum temperatures observed there have set new 
world records. In fact, the Antarctic winter season of 
1958 saw a continuous breaking and resetting of the 
vorld minimum temperature record. As it now stands, 
the -87.4° C. (—125.3° F.) recorded at Vostok on August 
25 is the world record, 

Now that this eventful winter is over, the still frag- 
mentary information from Antarctica can be analyzed 
ad summarized. However, a reservation must be made 
is to the reliability of sources of this information. We 
we compelled to resort to a practice usually avoided in 
wientific analysis; that is, to a utilization of newspaper 
ud radio reports. At this early date these are our only 
sources for some of the data which may be long in reach- 
ing publication i in technical journals. 

he a previous article [1] on the lowest temperatures 
reported on earth, we published the then latest record 
low air temperature for the world. It was —74.5° C. 
(-102.1° F.) observed at the U.S. Antarctic station at 


the South Pole (elev. 2,800 m.) on September 17, 1957. 
This minimum broke the previous record of —67.7° C. 
(—89.9° F.) observed at Oimekon in eastern Siberia on 
February 6, 1933. Carrying on from where we left off 
in the previous article, we show in table 1 a chronological 
list of the later world-record-breaking minima observed 
in Antarctica during the winter of 1958. 

The new record of —87.4° C. has exceeded a theoretical 
estimate of the lowest possible Antarctic temperature 
made by Shliakhov [2]. His estimate was based on the 
assumption that radiation balance becomes zero at the 
theoretical minimum temperature; i.e., the back radiation 
from the atmosphere equals the upcoming radiation from 
the surface. Using actual radiation measurements from 
Vostok 1 and theoretical values derived from soundings 
made at. Vostok 1 in May 1957 and at the South Pole in 
July 1957, and using the Shekhter radiation nomogram 
[3], Shliakhov computed the minimum temperature to 
vary from —78°+2° C. to —80°+1° C. He then per- 
formed another calculation based on greater station ele- 
vation (5 km. m.s.].) and using a mean lapse rate of 5° C. 
per 1,000 m. above the inversion layer as determined 
from radiosoundings made at the interior stations. He 
assumed the temperature of the snow surface to be 
—75° C., the strength of the inversion —30° C., the 
depth of the inversion layer 1 km., and also assumed that 


TARLE 1.—Temperature minima observed in Antarctica since 1957 


Station | Date Source 

Amundsen-Scott Station (South Pole)... —74.5 | —102.1 | Sept. 17, 1957_- E. C. Flowers and R. A. McCortrick, “World Record Low Temperature, South Pole, Septem- 

\, | ber 17, 1957,"" Monthly Weather Review, vol, 85, No. 11, Nov. 1957, p. 383. 

—108.4 | May 3, 1958__ Priroda Moscow, Aug. 1958, p. 82. 

Sovetskay: —110.2 | May 9, 10, 1958... Do. 

Vostok —113.8 | June 15, 1958__ Izvestiia, Moscow, July 2, 1958. 

Sovetskaya ok | —81.2 |) 114.2 | June 19; 1958... Do. 

Ostok _ —81.4 —114.5 | July 24, 1958____. .| Vodnyi Transport, Moscow, July 29, 1958. 

—83 —117.4 | July 25, 1958_____ Do. 

Ostok_ _ —85.5 | —121.9| Aug. 7, 1958....._| Vodnyi Tronepert Moscow, Aug. 19, 1958. 

ostok. _ | —85.8| —122.4] Aug. 8, 1958... M. J. Rub Rubin, “World Record Low Temperature,” Monthly Weather Review, vol. 86, No. 8, 

: ug. 1958, p. 308. 

Rubin, Record Low Temperature,” Monthly Weather Review, vol. 86, No. 8, 

958, p. 

Vostok. | —87.4 | —125.3 | Aug. 25, 1958_____- TASS, radio broadcast. The New York Times, Aug. 26, 1958. Izvestiia, Moscow, Aug. 

1958. E. 8. Rubinstein, “K Voprosu o Poliusakh Kholoda” (On the Problem of the Cold 
Poles), Meteorologiia i Gidrologiia, No. 12, 1958, pp. 28-30. E. M. Do’ hman, ““Opredekenie 
| Minimal’noi Temperatury v Antarktike” (Determination of the lowest temperature in 
| | Antarctica), Meteorologiia i Gidrologiia, No. 1, 1959, pp. 21-22. A. F. Treshnikov, Zakovannyi 
| Led (Icebound), Gosudarstvennoe Izdatel’ stvo Geografichesk oi Literatury, Moscow, 1959. 
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the air was saturated in all layers. Calculating the back 
radiation of the atmosphere (Z£’) from the Shekhter 
radiation nomogram and the corresponding equilibrium 
temperature (7',), he found: £’=0.100+0.005 cal./em.?/ 
min., and 7',=—85°+2° C. at 5 km. above sea level. 
Taking into account the above estimated lapse rate, 
Shliakhov came to the conclusion that, at the height of 
about 4 km., in Antarctica the lowest temperature could 
not be below —80°+2° C. At higher levels it could be 
lower by 0.5° per 100 m. 

The early low minima of May and June 1958 seemed to 
bear out this calculation, but the final record was 5.4° lower 
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than the predicted limit. In the following noie, Dy. 
Wexler attempts to explain why the theoretical minimyy 
temperature was too high. 
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NOTE ON LOWEST ANTARCTIC TEMPERATURE ESTIMATED 
BY SHLIAKHOV 


H. WEXLER 
U.S. Weather Bureau, Washington, D.C. 
(Manuscript received May 4, 1959] 


In the preceding note, Miss Stepanova referred to 
Shliakhov’s estimate of lowest Antarctic temperature [1] 
which has been proven too high by subsequent events. It 
appears that the basic error committed by Shliakhov is in 
the association of a zero radiation balance with a certain 
temperature—the theoretical minimum. This is not true 
within the likely range of atmospheric temperatures be- 
cause of the water-vapor absorption spectrum with its 
“window” near the maximum of the black-body radiation 
spectrum. As explained in an earlier paper [2], although 
the surface temperature may reach a state of “quasi- 
radiative equilibrium” wherein the back atmospheric 
radiation exactly equals the upcoming radiation from the 
surface, the surface temperature continues to fall for the 
following reason: All of the back radiation from the 
atmosphere is absorbed by the snow surface which acts as 
a black-body for infrared radiation but not all of the 
upcoming radiation is absorbed by the atmosphere which 
is not a black-body. Consequently there is a net loss of 
energy in the atmosphere which cools and then sends less 
back radiation to the surface, which then cools also, and 
On. 


Shliakhov essentially attempts to calculate this effect 
and comes to the conclusion that it could only change the 
temperature of the lowest 3 km. by 1° C. in 200 hours, or 
8.4 days, and neglects it in the calculation of the theoretical 
minimum temperature. But as the data in Stepanova’s 
table 1 show, the succession of minimum temperatures 
proceeds so slowly with time that a cooling rate of 1° C. 
per 8.4 days cannot be neglected. 

McCormick [3] has shown that in the absence of ad- 
vective and conductive processes the theoretical minimum 
could be —200° C. in 180 sunless days at the South Pole. 
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A NOTE ON VISIBILITY AT SACRAMENTO, CALIF. 


GEORGE C. HOLZWORTH 


=“ Weather Bureau Research Station, Berkeley, Calif. 
(Manuscript received January 13, 1959] 


ABSTRACT 


Daylight visibility observations at the Sacramento Municipal Airport were used to investigate the possibility 

’ that visibility in the area has been declining due to rising levels of air pollution. Those visibility observations 
made with naturally occurring fog (defined as relative humidity greater than 90 percent) and/or precipitation were 

eliminated from the data. Observations for July, a warm, dry month with relatively slight climatic differences 

from year to year, and November, a cool month that often exhibits a considerable annual variation of climate, were 

tabulated for the three 4-year periods, 1935-38, 1945-48, and 1953-56. The results show that for July the percent 

frequency of poor visibilities steadily increased, and of good visibilities steadily decreased over the last 20-odd 


years. 
lar. 
dispersion. 
annual climatic variations. 


Such a trend is consistent with the Sacramento population increase. 
This is believed to be due to variations in the occurrence of meteorological factors important to atmospheric 
It is likely that such irregularities are characteristic of other months and areas where there are large 


For November the trend was irregu- 


From concurrent surface wind observations, the preferred wind conditions for various visibility ranges were 


determined. 


These data indicate that poor visibilities at the airport, on the southern edge of the city, occurred 


most frequently with light winds from over the nearby Sacramento urban area. Good visibilities were favored 


by moderate winds from rural directions. 
for loose materials to be picked up from the ground. 


1. INTRODUCTION 


Residents of Sacramento and other cities in California’s 
Central Valley, as well as aircraft pilots in the vicinity, 
have stated that over the years they have observed a gen- 
eral decrease in visibility. These individuals have at- 
tributed such a trend to increasing levels of air pollution. 
Because air pollution measurements have not been made 
in the Central Valley until recently, studies of trends in 
the concentrations of specific pollutants are not possible. 
However, the solid and liquid particulate matter sus- 
pended in the atmosphere, which represents air pollution, 
reduces visibility [1], and therefore a measure of visibility 
may be interpreted as a measure of air pollution. Hence, 
changes in visibility not attributable to naturally occurring 
aerosols may be indicative of manmade trends in air pol- 
lution. This study was undertaken in order to assess the 
possibility that there has been a general deterioration of 
visibility in the Sacramento area due to human activity. 


2. DATA AND TREATMENT 


From the hourly weather observations made by the U.S. 
Weather Bureau at Sacramento Municipal Airport, the 
July and November data for the three 4-year periods 


Visibilities were reduced when the wind speed became strong enough 


1935-38, 1945-48, and 1953-56 were selected and punched 
on cards for machine tabulation. Only daytime observa- 
tions, 0730-1730 psr,' were considered since visibility at 
night depends upon different factors and the observations 
are generally less reliable [1]. All observations, 1935-56, 
were made from the same location. 

Prior to 1943 visibility to be reported by Weather 
Bureau observers was defined as “the greatest cistance 
toward the horizon that known objects can be identified.” 
This was later revised to “the lowest visibility over the 
half of the horizon with the greatest visibility.” Sum- 
maries for this study do not suggest that this change in 
definition had any noticeable effect on Sacramento 
observations. 

July and November visibility trend analyses were con- 
sidered with respect to the three 4-year periods. The data 
were lumped into 4-year groups in an attempt to smooth 
out year-to-year variations in meteorological conditions 
such as factors affecting the dispersion of pollutants. 

The natural aerosols which limit visibility are fog and 
precipitation. Precipitation may occur at various relative 


1 Actual observation times are about 28 minutes past the hour for 194- 
56 and about 41 minutes past the hour for 1935-38. At Sacramento in 
November the latest sunrise is about 0700 pst and the earliest sunset i8 
near 1700 Pst. 
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MILES 7-10 MILES II-IS MILES 20-30 MILES >30MILES 


Figure 1.—July visibilities at Sacramento Municipal Airport, 0730— 
1730 pst. Percent frequencies when relative humidity less than 
#1 percent and no precipitation occurring. 


humidities, but natural fog is unlikely below about 91 
percent. Those visibility observations with precipitation 
and/or relative humidity greater than 90 percent were 
tabulated separately, and remaining visibility observations 
were originally broken down according to various other 
relative humidity ranges. Although the visibilities tended 
to be lower with higher humidities, the trends between 
tyear periods for various humidity ranges were all 
similar. For this reason, and because it is desirable to 
deal with large samples, all visibilities recorded with 
humidities below 91 percent (and no precipitation) were 
considered together. In 96.7 percent of the July visibility 
cases, humidities were less than 71 percent. In November, 
humidities were below 71 percent in 60.2 percent of the 
cases and below 91 percent in 83.9 percent of the cases. 


3. RESULTS 


In the following discussions visibility observations with 
precipitation and/or relative humidity greater than 90 
percent are excluded. This provision is intended to elimi- 
nate visibility restrictions that might be attributed to 
naturally occurring fog and precipitation particles. 


JULY 


July was chosen as typical of the warm, arid summers 
when there are seldom marked variations in Central Val- 
ley weather. The percent frequencies of visibilities in 
given ranges for each of the three 4-year periods are shown 
infigure 1, This bar diagram indicates that the frequency 
of visibility reports in the ranges 3-6 and 7-10 miles 
steadily increased over the years studied. The percentage 
frequency of the combined 3—10-mile range increased by 
4 factor of almost 5 from 1935-38 to 1953-56, apparently 


at the expense of the 20-30-mile range whose frequency 
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11935 -38 1945-48 


1953 - 
PERCENT (953-56 


Figure 2.—November visibilities at Sacramento Municipal Airport, 
0730-1730 pst. Percent frequencies when relative humidity less 
than 91 percent and no precipitation occurring. 


steadily declined from 62.8 percent in 1935-38 to 44.5 per- 
cent in 1953-56. The frequency of visibilities of 11-15 
miles changed only slightly and inconsistently with time. 
In the range exceeding 30 miles, there was no regular 
trend. There were no July visibilities of 0-214 miles. 
Considering all visibilities, a deteriorating trend is clearly 
indicated for July. 

The most striking visibility trends were in the 3-10- 
and 20-30-mile ranges. By the chi-square test [2], dif- 
ferences between the three 4-year periods for both these 
visibility ranges are significant at the 1-percent level. It 
should be mentioned that such tests of significance assume 
independence of samples, and hourly weather observations 
do not fully satisfy this assumption. However, each of 
the 4-year groups was considered to more clearly approach 
the property of independence. 


NOVEMBER 


Climatologically, November is a relatively cool month 
and was chosen to represent early winter conditions as 
opposed to summer. The weather is often quite variable 
within this month and considerable differences may occur 
from year to year. Some Novembers are rainy and foggy 
while others are relatively dry and fog free. 

The percentage frequencies of visibilities in given ranges 
for each of the three 4-year periods, as shown in figure 2, 
indicate that November visibilities followed an irregular 
trend. As a whole, poor visibilities (0-10 miles) were 
least frequent, and good visibilities (20 miles or greater) 
were most frequent in the middle period, 1945-48. In the 
later period, 1953-56, the picture was reversed with 0-10- 


mediate between those of the other two periods. The fre- 
quency of visibilities of 11-15 miles differed bet ween peri- 
ods by only 0.2 percent. Considering all visibilities for 
November, they may be described as having been inter- 
mediate in 1935-88, good in 1945-48, and poor in 1953-56. 
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This irregular variation stands in sharp contrast to the 
steady trend shown for July. As mentioned above, data 
for 4-year periods were lumped together in the hope that 
important weather variations would be smoothed out. 
This result is more probable for July than for November 
when conditions may be quite variable. November aver- 
age 4-year precipitation amounts were 1.09, 1.41, and 1.59 
inches for the earliest to latest periods, respectively ; for 
any one November in these 12 years, precipitation ranged 
from 0.03 to 3.35 inches. These are cited to indicate to 
some extent the degree of weather variability in these 12 
Novembers. It is thought that such variations in impor- 
tant meteorological factors, affecting both dispersion and 
emission of visibility-reducing pollutants, have resulted 
in the irregular November trend shown in figure 2. 
With respect to pollutant sources, the burning of crop 
residues after harvesting should be mentioned. It is 
thought that this practice, which has been common in the 
rich agricultural area of Sacramento and other parts of 
the Central Valley, is partly responsible for the irregular 
isibiliti Smoke from what agri- 
cultural burning occurs in July is likely to be more rapidly 
dispersed by the usual summer conditions of windiness 
and instability. However, because the desired facts on 
agricultural burning are nonexistent, such remarks must 
be regarded as speculation. 


4. FURTHER CONSIDERATION 


It is generally recognized that as urban areas expand, 
pollutant emissions from such sources also increase (in 
the absence of emission controls). It is logical then that, 
other things being the same, visibility should show a gen- 
eral decline with increasing urbanization. The Sacra- 
mento Municipal Airport is on the southern edge of the 
city, about 4144 miles from the downtown area. Because 
Sacramento is the only large city in the county and because 
much of the total metropolitan population lies beyond the 
official city limits, Sacramento County population figures 
have been used. From decennial population data [3] and 
others [4], the following population estimates were made 
for Sacramento County for 1 year within each of the 


three periods studied : 


associated with a rapidly increasing population. In No- 


vember this same trend held true for 1945-48 to 1953-56, 
but improving visibilities from 1935-38 to 1945-48 were 
also associated with a population increase. As pointed 
out, this latter irregularity is believed to have been due to 
variations of important meteorological factors. 

That the Sacramento metropolitan area was indeed a 
source of visibility-reducing pollutants may be shown from 
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Fieure 3.—Schematic map of Sacramento, Calif. 


a consideration of surface winds. For each visibility ob- 
servation the concurrent wind observation was tabulated 
in wind rose form for each visibility range. Wind diree- 
tions were divided into two groups, northwest clockwise 
through east-southeast (NW-ESE) and southeast through 
west-northwest (SE-WNW). As shown in figure 3, most 
of the Sacramento metropolitan population lies in the 
sector NW-ESE from the airport. The other sector, SE- 
WNW, may be described as relatively rural. The rela- 
considered first. This was done by summarizing the per- 
cent frequency of occurrence of wind in each of the two 
direction sectors, NW-ESE and SE-WNV, by 5-m.p.h. 
bilities. The ratio of the former to the latter was then 
determined and these values for July and November are 
given in table 1. Each number in this table is the ratio 
of percent frequency of given wind conditions for low 
visibility to the percent frequency of such wind conditions 

isibiliti Thus, the higher the ratio, the greater 


from NW-ESE at 1-5 m.p.h. It should be made clear 
that these wind conditions are not necessarily those with 
which the visibility was most often 0-10 miles. The pref 
erence is relative to what is expected from a consideration 


2 Using tables B and F of Local Climatological Data Supplement, it 8 
possible to compute the ratios of “percent of time the wind is from 
given directions for visibilities less than 7 miles with smoke or haze” to 
“percent of time the wind is from corresponding directions for all 
visibilities.”’ 


&§ 
“ene 
the “preference” of iow visibilities for such wind condi- 
tions.2- For example, table 1 shows that in July the most 
&§ preferred wind condition for 0-10-mile visibilities was 
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Taste |.—Ratio of “percent of winds for visibility 0-10 miles’’ to 
“percent of winds for all visibilities.’’ Data for Sacramento 
Municipal Airport when relative humidity less than 91 percent and 
no procipttation ; 0730-1730 pst, 1935-38, 1945-48, 1953-56 
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TaBLeE 2.—Ratio of ‘“‘percent of winds for visibility 20-80 miles’’ to 
“percent of winds for all visibilities’’. Data for Sacramento Munic- 
ipal Airport when relative humidity less than 91 percent arid no 
precipitation; 0730-1730 pst, 1935-38, 1945-48, 1953-56 


Wind direction Wind direction 
Winds Wind s 
m.p.n. 
NW-ESE |SE-WNW| Calm All NW-ESE |SE-WNW]| Calm All 
JULY JULY 
(No visibility reported in 0-242 miles range) 

NOVEMBER 

NOVEMBER 
0. 34 0. 34 

1. 36 1.16 1.33 2 30 1. 86 
-9 . 89 


Wind speeds exceeding 15 m.p.h. were infrequent, and 
so the ratios for such speeds had to be computed from 
small percentages. However, the computed ratios appear 
to exhibit reasonable trends as wind conditions vary. A 
likely exception is the value of 2.23 in table 1, July, which 
probably would be lower with more cases of such winds. 
This fault might be corrected by combining some of the 
higher wind speeds, but in so doing one of the features 
of the table would be lost. It is a reasonable hypothesis 
that, considering nearby community sources of visibility- 
reducing aerosols, improving visibility should be antici- 
pated with increasing wind speed (more favorable dilu- 
tion) up to the point where dust and other loose materials 
are picked up by the wind. This is borne out in table 1 
where, considering all wind directions, there is a decrease 
in the ratios with increasing speeds up to the point where 
loose material is picked up. Here the ratios begin to 
increase. In July when precipitation averaged less than 
a trace and the ground was dry, this reversal began at 
16-20 m.p.h. In November when precipitation averaged 
136 inches, the ratios did not begin to increase until 
speeds exceeded 20 m.p.h. 

From table 1 it is clear that wind directions occurring 
area. This was especially true for July when for every 
wind speed the ratios were greater for the sector NW- 
ESE than for the SE-WNW sector. For November, con- 
sidering all wind speeds, the most preferred direction with 
low visibilities was also from the metropolitan area. 
However, the preference was not nearly so strong as in 
July. This may have been due in part to agricultural 


burning which provides pollutant sources to the south 
as well as to the north. 

Table 2 is similar to table 1 except it is for good visi- 
bilities, 20-30 miles. This table indicates that for all 
directions there was an increase in the ratios, preference 
point where loose materials are picked up, and then the 
trend reversed. Considering all wind speeds, the most 
preferred direction for good visibility was from the rural 
sector, SE-WNW. 

Hygroscopic nuclei may be active in the humidity 
range 71-90 percent. If such relative humidities are re- 
lated to wind directions, then the concurrent visibilities 
might be lower and tables 1 and 2 would be invalid. July 
daytime humidities above 70 percent are very rare. In 
November, humidities of 71-90 occurred 23.7 percent of 
the time. The November ratios in tables 1 and 2 were 
recomputed, using only data for relative humidities be- 
low 7i percent (and no precipitation). The recomputed 
ratios are almost identical with the original. It is con- 
cluded that any relationship between humidities of 71-90 
percent and wind directions does not alter the ratios of 
tables 1 and 2. 


5. CONCLUSIONS 


have been deteriorating over the last 20-odd years. This 
is most apparent in the 0-10-mile visibility range, which 
showed a steady increase in percent frequency, and in 
the 20-30-mile range, which showed a steady decrease. 
from intermediate to good to poor, from the earliest to 
the latest periods, respectively. This may have been due 
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largely to variations in the occurrence of meteorological 
factors, especially those important to atmospheric dis- 
persion of manmade pollutants. It is likely that such 
irregularities are characteristic of other seasons and areas 
where there are large variations in the weather. 

From a consideration of wind data it is clear that poor 
visibilities were associated with wind directions from over 
the urban area and with light winds. Wind direction and 
speed factors were interdependent. Good visibilities were 
most likely with winds from rural directions and at 
reduced when winds became strong enough to pick up 
dust and other loose materials. 
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THE WEATHER AND CIRCULATION OF APRIL 1959 
Including the Role Played by an Index Cycle 


JAMES F. ANDREWS 
Extended Forecast Section, U.S. Weather Bureau, Washington, D.C. 


1. INTRODUCTION 


April 1959 was characterized by a lack of severe local 
storminess throughout most of the United States. A sim- 
ilar condition prevailed in April of 1958 [1], in sharp 
contrast to the stormy April of 1957 [2] when a record 
number of tornadoes was observed. From the northern 
Rocky Mountain States southeastward to the Gulf States, 
temperatures averaged below the seasonal normals, while 
west of the Rockies and in the northeastern third of the 
Nation the weather was generally mild. The long period 
of unusual warmth continued along the California coast, 
with new April temperature records being established. 
Precipitation continued deficient in the northern Great 
Plains and upper Mississippi Valley, with some areas re- 
porting the least amount of record for the period January 
to April. Subnormal amounts of precipitation also fell 
in parts of the Southwest, thus continuing the dryness of 
the winter season. bg 

Despite the fact that few monthly records were estab- 
lished during April, pronounced changes in weather and 
circulation did occur. These changes can be related to a 
strong index cycle. 


2. MONTHLY MEAN CIRCULATION 


The average circulation pattern for April 1959 (fig. 1) 
was one of relatively fast westerly flow with wave systems 
of small amplitude around the entire Northern Hemi- 
sphere. Over the Western Hemisphere the zonal index, a 
measure of the net geostrophic west wind between latitudes 
5° N. and 35° N., averaged 0.8 m.p.s. above the April nor- 
mal. This continues the period of above normal values of 
this index that began in February 1959 [3]. At higher lati- 
tudes the polar westerlies (55° N.-70° N.) were 0.8 m.p.s. 
below normal, mostly as a result of blocking over eastern 
Canada. The subtropical westerlies (20° N.-35° N.) also 
aweraged slightly below their normal value. 

Only one full-latitude trough was observed during 
April, from the Bering Sea southward to just west of the 
Hawaiian Islands (fig. 1). Nearly-full-latitude troughs 
Were observed along the eastern Asiatic coast and over 
eastern North America. The latter trough, along with a 


strong ridge over the Pacific coast, controlled much of the 
weather in the United States. Blocking in Scandinavia 
was associated with a middle- and low-latitude trough 
extending from eastern Europe to North Africa. The 
latter two features of the planetary circulation resulted in 
confluence over central Asia, where a narrow band of above 
normal wind speeds prevailed at 700 mb. 

In the Pacific the zone of maximum wind speed at T00 
mb. was found just east of the Japanese Islands (fig. 2A), 
somewhat farther west than normal. This position can 
be related to the westward displacement of the normal 
trough near the Asiatic coast [4] and to the appearance 
of a mean ridge in the western Pacific (fig. 1). 

Across North America the major trough-ridge system 
was also displaced west of its position on the normal map 
[4], though wind speeds at 700 mb. over the United States 
did not depart significantly from their normal values 
(fig. 2B). As a result the observed jet axes at 700 mb. 
bore a strong resemblance to their normal counterparts 
(not shown). Fast westerly flow. prevailed across the 
Atlantic, where wind speeds were above normal in a nar- 
row band centered around the 45th parallel (fig. 2). 
Greatest departures were observed in the eastern Atlantic, 
where speeds were as much as 6 m.p.s. above normal. 
Here also was found the center of greatest height 
anomaly over the Western Hemisphere, some 260 ft. 
below normal (fig. 1). Mean sea level pressures in this 
area were as much as 8 mb. below normal. A well-defined 
jet axis existed across the Atlantic (fig. 2A). This cur- 
rent split into two branches as it approached Europe, one 
branch passing northward through Scandinavia, the other 
southward across North Africa. Here again the con- 
figuration of the primary jet axes strongly resembled the 
normal. 

The greatest anomaly of 700-mb. height, —310 ft., was 
associated with an extensive center of action over northern 
Asia (fig. 1). Very cold air, as much as 10° C. below 
normal in the layer from 1,000 mb. to 700 mb., accom- 
panied this center. Because of this extreme cold, sea 
level pressures, on the average, were at most only 4 mb. 
below normal. 
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Fieure 1.—Mean 700-mb. height contours (solid) and departures from normal (dotted) (both in tens of feet) for April 1959. Prin- 
cipal circulation features influencing weather in the United States were stronger than normal ridge near the west coast and 


trough in eastern North America. 


3. AVERAGE UNITED STATES WEATHER 
RELATED TO THE MEAN CIRCULATION 


Charts depicting the departure of average surface tem- 
perature from normal and the total precipitation for 
April are shown in figure 3. In addition, figure 4 gives 
the percentage of normal precipitation for April and 
for the period December 1958—April 1959. 

A warm dry regime prevailed quite generally west of 
the Continental Divide, with temperature departures of 


TABLE 1.—Selected California cities which reported unusual warmth 
during April 1959 


Mean tem-| Anomaly} Years of 
City rs a (° F.) record Remarks 
ngeles (city office) 

Los Angeles 65.4) +6.2 28 ||Warmest April of 
Sa 64.7 +4.2 87 record. 
San Francisco (airport) --__-__- 59.8 +5.7 31 
60.1 +4.5 31 of 
SS 64.3 +5.1 82 ye warmest April 
51.8 +6.7 54 record. 
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Fieure 2.—(A) Mean 700-mb. isotachs and (B) departure from 
normal wind speed (both in meters per second) for April 1959. 
Solid arrows in (A) indicate primary axes of westerly jet, with 
speeds greater than 12 m.p.s. stippled. Jet maxima were located 
hear the Asiatic coast and over the Atlantic Ocean, while wind 
speeds were close to normal over most of the United States. 


as much as 6° F. in portions of Arizona and California 
(fig. 3A). In the latter State the warmth was of a 
record-breaking nature, thus continuing the long period 
of unusual warmth which began during winter 1957- 
198. Table 1 shows a representative list of cities in 
California where heat records were approached or estab- 
lished. Precipitation was markedly deficient throughout 
much of the Far West, with approximately 25 percent 


Fieurn 3.—(A) Departure of average surface temperature (° F.) 
from normal and (B) total precipitation (inches) for April 
1959. Features were the record warmth along the California 
coast and continued precipitation deficiency in the northern 
Plains. (From Weekly Weather and Crop Bulletin, National 
Summary, vol. XLVI, No. 18, May 4, 1959.) 


of normal in California, Nevada, southern Idaho, and 
western Arizona (fig. 4A). Since December 1958 a vast 
area from central Texas to southern California has re- 
ceived less than half its normal amount of precipitation 
(fig. 4B). Amarillo and Wichita Falls, Tex., have been 
drier than normal for 9 consecutive months. Abundant 
sunshine accompanied the warm, dry weather, with Ely, 
Nev., and Sacramento, Calif., experiencing their sunniest 


April on record. 


Anticyclonic conditions associated with a stronger than 
normal upper-level ridge near the west coast (fig. 1) were 
largely responsible for this warm, dry regime. Moreover, 
northerly anomalous flow tended to inhibit precipitation. 
In California stronger than normal northeasterly flow 
resulted in foehn warming, thus contributing further to 
the unusual warmth along the coast. 

Temperatures averaged below normal in the Great 
Plains, the lower Mississippi Valley, and along the Gulf 
coast (fig. 3A). Greatest departures were observed in 
southeastern Texas and northwestern Louisiana, where 
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Figure 4+.—Percentage of normal precipitation for (A) April 1959 
and (B) December 1958—April 1959. Note prolonged period of 
dryness in Southwest and northern Plains. (From Weekly 
Weather and Crop Bulletin National Summary, vol. XLVI, No. 
18, May 4, 1959.) 


temperatures were as much as 4° F. below normal. Most 
of this coolness can be related to polar air masses associated 
with anticyclones steered southward by stronger than 
normal northerly flow (fig. 1) along a_ well-defined 
primary track (fig. 5B). 

Precipitation in the Great Plains and Mississippi Valley 
was generally well below normal (fig. 4A). A large area 
in the Dakotas and Minnesota received less than a half 
inch of moisture during April (fig. 3B). In that area 
the total precipitation for the 5-month period December 
1958 to April 1959 has been much below normal (fig. 4B). 
At Duluth and St. Cloud, Minn., with 1.34 and 1.22 inches 
for the period January to April, it was the driest such 
period in their respective 88 and 54 years of record. At 
Huron, S. Dak., this same period was the second driest 
in 77 years. 

Most of this precipitation deficiency in the Plains States 
and Mississippi Valley appears to have been a function of 
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Ficure 5.—Frequency of (A) cyclone passages and (B) anti- 
cyclone passages (within equal area boxes of 66,000 n. mi.’ [6]) 
during April 1959. Primary tracks are indicated by solid ar 
rows. Cyclone tracks in the United States were rather well 
related to the precipitation pattern, while the major anticyclone 
track was associated with cool weather in mid-United States. 


stronger than normal northwesterly flow to the rear of 
the mean trough in eastern North America (fig. 1). In 
addition, the primary storm tracks passed both north and 
south of the northern Great Plains and upper Mississippi 
Valley (fig. 5A). Along the southern track a narrow band 
of heavier precipitation extended from southern Montana 
through Nebraska to the lower Lakes Region (figs. 3B and 
4A). In Wyoming and Colorado, precipitation was largely 
in the form of snow which occurred during two periods 
and brought the month’s total snowfall generally well 
above average. Heavy precipitation also occurred along 
the Texas coast, mostly as the result of overrunning of coo 
polar air masses in southwesterly flow aloft. 
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Fieurr 6.—Time variation of speed of 700-mb. westerlies averaged 
over the Western Hemisphere between latitudes 35° N. and 
55° N. Solid line connects 5-day mean index values (plotted at 
middle of period and computed thrice weekly), while dashed 
line gives the corresponding normal. Index cycle was a 
prominent feature of the circulation. 


Above-normal temperatures prevailed quite generally 
east of a line from the central Dakotas to southern Georgia, 
with greatest departures in the Middle Atlantic States 
(fig. 3A). This unseasonable warmth was associated with 
above normal heights and southerly anomalous flow at 
700 mb. (fig. 1) and at sea level (chart XI in [5]). The 
heaviest precipitation in the United States occurred rather 
generally throughout the eastern half of the Nation (fig. 
3B) and was related primarily to the mean trough (fig. 1). 
At sea level two principal storm tracks, one across the 
Great Lakes and the other along the Atlantic coast (fig. 
5A), contributed to this precipitation. Lakeland, Fla., 
reported an alltime April precipitation record of 8.48 
inches, with the January—April total of 28.00 inches also a 
record for the period. A relatively short distance away, 
Key West, Fla., with a long station record, received no 
measurable precipitation, also an April record. This ex- 
treme difference between Lakeland and Key West cannot 
readily be explained by large-scale, time-averaged charts. 


4. WEEK-TO-WEEK VARIABILITY 


Figure 6 shows the time variation of the temperate- 
latitude 5-day mean zonal index for April 1959. The com- 
plete index cycle portrayed here was accompanied by 
pronounced changes in weather and circulation. These 
changes will be described, using a series of 5-day mean 
700-mb. charts centered 1 week apart, along with accom- 
panying charts of weekly temperature departure from 
normal and total precipitation (figs. 7 to 11). 


WEEK ENDING APRIL 5 


The month began with the westerlies at their greatest 
‘trength, some 4 m.p.s. above normal (fig. 6), with fast 
flow predominating over nearly the entire hemisphere 
(fig. 7A). The only exception to this strong cireumpolar 
whirl was in Europe, where a strong blocking ridge, ac- 
companied by a Low in the Mediterranean, resulted in a 


4, Departure of Average Tex.sperature from Normal for the Week Ending 


Total Precipitation, Inches, for the Week Ending 
a Midnight, | st, April 5, 1959 


Based on preiiminary teregrapnic reports 


Fieure 7.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for 
March 31—April 4, 1959. (B) Departure of average temperature 
from normal (° F.), and (C) total precipitation (inches) for 
week ending April 5, 1959. (B and C from Weckly Weather 
and Crop Bulletin, National Summary, vol. XLVI, No. 14, April 
6, 1959.) 
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Y Departure of Average Temperaiure from Normal for the Week Ending 
A 


Based on preliminary telegraphic reports 


Fieure 8.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for 
April 7-11, 1959. (B) Departure of average temperature from 
normal (° F.), and (C) total precipitation (inches) for week 
ending April 12, 1959. (B and C from Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVI, No. 15, April 13, 
1959. ) 
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more meridional flow pattern. The strong resemblance 
between the 5-day mean 700-mb. chart for the first week 
(fig. TA) and the mean 700-mb, chart for the month (fig, 
1) is noteworthy. 

Unseasonably warm weather accompanied this high 
index circulation over nearly the entire United States 
(fig. 7B). Mild Pacific air masses, warmed by foehn ae- 
tion in the Rocky Mountains, dominated the country and 
resulted in record or near-record daily maximum tem- 
peratures throughout much of the northern and central 
Rocky Mountain States. For instance, Salt Lake City, 
Utah, reported a record 82° F. on April 5, while Pueblo, 
Colo., with 83° F. on the 2d, equaled its previous April 
record. 

Precipitation during the first week was confined mostly 
to the extreme Northwest and the area from the middle 
and lower Mississippi Valley eastward (fig. 7C). Most 
of the latter precipitation was associated with the deep 
trough in the Ohio Valley (fig. 7A) and fell as a deepen- 
ing storm moved northeastward from the lower Great 
Plains (chart X in [5]). Tornadoes on the 2d caused 
some loss of life and extensive damage in northeastern 


Texas. 
WEEK ENDING APRIL 12 


Gradual amplification of the circulation pattern oe- 
cured during the second week as the westerlies diminished 
to more normal values (fig. 6). In the mid-Pacifie a 
stronger than normal ridge appeared, while farther east 
the mean trough of the previous week deepened while 
remaining nearly stationary (figs. 7A and 8A). Greatest 
changes occurred in North America, where strong ridg- 
ing and a large area of positive height anomaly developed 
along the Pacific Coast (fig. 8A). At the same time the 
trough in eastern North America sheared, the higher 
latitude portion moving eastward, while the low-latitude 
portion retrograded and redeveloped in the Southwest. 
Blocking began to appear in the North Atlantic as an 
intense cyclonic center of action shifted from the eastern 
Greenland coast to the North Sea, thus replacing the 
strong ridge of the week before (figs. 7A and 8A). This 
rapid change in circulation was accompanied by an 
equally sharp change in weather throughout western 
Europe—from sunny and pleasant to cold and stormy. 
Blocking also appeared in the polar regions in the form 
of a high pressure area with 700-mb. heights as much as 
630 feet above normal. 

In the United States the change in circulation broughit 
a rapid return to winter weather over much of the Nation. 
Stronger than normal northerly flow (fig. 8A) deployed 
very cold polar air into the United States, where tem 
peratures for the week averaged well below norma! from 
the Rocky Mountain States to the Ohio Valley—as much 
as 15° F. below normal at Denver, Colo. (fig. 8B). The 
cold was ushered in on the heels of a deep Alberta-type 
disturbance whose strong winds created severe dust storms 
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in the Dakotas, Minnesota, and Wisconsin on the 6th and 
7th. Many new daily minimum temperature records were 
established by the abnormal celd in the northern and 
central Rockies, where only a few days before record 
warmtl had prevailed. Some of these records were (in 
oF.): 12° at Helena, Mont., on the 9th; 7° at Denver, 
Colo., on the 10th and 12th; and 6° at Pueblo, Colo., on 
the 10th. In addition, a temperature of —3°F. at 
Cheyenne, Wyo., on the 12th was the coldest ever observed 
there so late in the season. The Far West and the East, 
where 700-mb, heights were above normal (fig. 8A), con- 
tinued warm (fig. 8B). The period of greatest warmth 
in the East, from the 7th to 9th, saw numerous daily maxi- 
mum temperature records established from Kentucky to 
southern New England. Baltimore, Md., reported the 
highest temperature, 92° F., on the 9th. 

Precipitation was mostly moderate from Texas north- 
eastward, with little or none over the Far West and 
northern Great Plains (fig. 8C). Snows were frequent 
over the Rocky Mountains and extended from the Texas 
Panhandle to the western Great Lakes. Over a foot of 
snow fell at the higher elevations in Wyoming and Colo- 
rado. Rains were rather persistent in Texas, with heaviest 
amounts along the upper coast (fig. 8C). The gradual 
eastward push of the cold air was accompanied by up to 
6 inches of snow in parts of West Virginia, Maryland, and 
Pennsylvania on the 12th. The 6-inch fall at Charleston, 
W. Va., was the greatest snowfall ever recorded there in 
April. 

This precipitation pattern was the result mostly of 
overrunning of cold polar air by warmer, more moist, 
tropical air in southwest flow aloft (fig. 8A).~ At sea 
level, high pressure dominated the North American conti- 
nent, with easterly winds prevailing across most of the 
United States. The polar front lay in a nearly stationary 
position from the Middle Atlantic States to the western 
Gulf of Mexico for nearly the entire week. 


WEEK ENDING APRIL 19 


The westerlies continued to diminish, with their aver- 
age strength below normal during this period (fig. 6). 
At 700 mb, the planetary circulation pattern assumed 
many of the characteristics commonly associated with a 
low zonal index. These include a relatively short wave 
spacing between trough-ridge systems, amplification, and 
increased cellular nature of the pattern (fig. 9A). The 
latter feature, which also characterizes blocking, is further 
illustrated by the 700-mb. height anomaly field. Note 
(fig. 9.) the preponderance of positive anomaly at higher 
latitudes and negative anomaly at lower latitudes over 
much of the Western Hemisphere. 

The trough near the Asiatic coast moved eastward 
during the week, while retrogression of the ridge in 
western North America to the Gulf of Alaska cut off the 
eastern Pacific trough and limited it to a deep system 
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Ficure 9.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for 
April 14-18, 1959. (B) Departure of average temperature from 
normal (° F.), and (C) total precipitation (inches) for week 
ending April 19, 1959. (B and C from Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVI, No. 16, April 20, 
1959.) 
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at. lower latitudes (fig. 9A). As a consequence, retrogres- 
sion also occurred in North America, where a nearly-full- 
latitude trough appeared in the West, while ridging over 
the East replaced a strong confluence zone. As blocking 
spread westward and southward from its position over 
polar regions the previous week, it became entrenched 
from the Davis Strait to Europe (figs. 8A and 9A) and 
suppressed the fast westerly belt over the Atlantic and 
southern Europe. 

Temperatures in the United States remained unseason- 
ably cool during the third week over all but the Far South- 
west and Northeast (fig. 9B). The cool weather in the 
West was related to northwesterly flow aloft along the 
Pacific coast and below normal heights (fig. 9A). In 
the Gulf States, where heights were close to normal, the 
below-normal temperatures accompanied a very cold 
anticyclone which moved from Texas slowly eastward 
(chart IX in [5]). The warmth in the Northeast can be 
attributed to the stronger than normal mean ridge (fig. 
9A). 

Precipitation was generally heavy from Wyoming to 
Iowa and from Oklahoma to North Carolina (fig. 9C). 
Elsewhere amounts were mostly light to moderate, except 
in California and portions of the western Plateau where 
no measurable amounts fell. Snow fell once again from 
the Central Rockies to the lower Great. Lakes late in the 
week, with a total fall of 15 inches at Lander, Wyo. A 
series of waves developing along a cold front as it moved 
slowly eastward and southeastward across the Nation was 
responsible for much of the precipitation during the week. 


WEEK ENDING APRIL 26 


The zonal index reached a minimum value of 6.5 m.p.s. 
between the third and fourth weeks, and then, although 
remaining below normal, began to recover (fig. 6). The 
principal seat of blocking spread from the Davis Strait 
to the Pacific with characteristic pronounced meridional 
flow pattern, strong ridges, and cut-off lower latitude Lows 
(fig. 10A). As blocking relaxed in Canada, the westerlies 
increased across North America, the trough in the south- 
western United States moved rapidly eastward across 
the country and amalgamated with the trough in the 
western Atlantic, while the ridge off the Pacific coast also 
moved eastward to the Rocky Mountains (figs. 9A and 
10A). 

Subnormal temperatures predominated over most of the 
United States (fig. 10B) as stronger than normal north- 
westerly flow in western Canada (fig. 10A) continued to 
deploy cold Canadian air masses into the country. How- 
ever, eastward motion of the trough-ridge system in the 
United States was related to a warming trend which grad- 
ually spread eastward. This is readily seen by comparing 
figure 9B with figure 10B. Note the growth of the area 
of above normal temperatures in the Far West, associated 
with ridging and above normal heights, and the diminu- 
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Figure 10.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for 
April 21-25, 1959. (B) Departure of average temperature from 
normal (° F.), and (C) total precipitation (inches) for week 
ending April 26, 1959. (B and © from Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVI, No. 17, April 2%, 
1959.) 
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“ht Fieure 11—(A) 5-day mean 700-mb. contours (solid) and height 
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tion of the area of above normal temperature in the North- 
east, related to the decrease in 700-mb. height (figs. 9A 
and 10A). At Roswell, N. Mex., the temperature reached 
97° F. on the 25th, an April record. 


The precipitation pattern also tended to move eastward 


from the preceding week (figs. 9C and 10C). Heaviest 


a 


mounts fell in the eastern Gulf States and portions of 


the Atlantic Coastal States. Most of this occurred early 
in the week in connection with a cold front and cyclonically 


Cc 


urved contours at 700 mb. (fig. 10A). As the west coast 


ridge moved inland and the eastern Pacific trough ap- 
proached the coast, precipitation broke out in California, 
southern Oregon, and western Nevada. Rain on the 25th 
at Burbank, Calif., ended a period of 62 consecutive days 
without measurable precipitation. 


WEEK ENDING MAY 3 


A peak value of 9.6 m.p.s. was reached by the zonal index 


between the fourth and fifth weeks, thus completing the 
cycle which began at the beginning of April (fig.6). This 
increase in the mid-tropospheric zonal circulation was 
manifested at the 700-mb. level by eastward motion of 
nearly all the mean trough and ridge systems (figs. LOA 
and 11A). 


The warming trend in evidence in the West the previous 


week overspread almost the entire country, with only the 
Northwest and extreme Northeast experiencing below nor- 
mal temperatures (fig. 11B). Greatest departures were 
reported in the middle of the Nation, where a large anom- 


a 
c 


ly area greater than +9° F. was observed. Heavy pre- 
ipitation, the result of strong onshore cyclonic flow, re- 


turned to the Pacific Northwest (fig. 11C), while heavy 
amounts were also reported from Iowa to southern New 
England as a by-product of the trough over the Great 
Lakes (fig. 11A). 


r 


April came to a close with the westerlies decreasing 
apidly, and as May began they were once again well 


below normal (fig. 6). 
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Suggestions For Authors 


Articles are accepted for the Monthly Weather Review 
with the understanding that they have not been published 
or accepted for publication elsewhere. 


Two copies of the manuscript should be submitted. All 
copy, including footnotes, references, tables, and legends 
for figures, should be double spaced with margins of at 
least 1 inch on sides, top, and bottom. Some inked cor- 
rections are acceptable, but pages with major changes 
should be retyped. The style of capitalization, abbrevia- 
tions, etc., used in the eview is governed by the rules set 
down in the Government Printing Office Style Manual. 


Tables should be typed, each on a separate page, with 
a title provided. They should be numbered consecutively 
in arabic numerals. 


In equations, conventional symbols in accordance with 
the American Standards Association Letter Symbols for 
Meteorology should be used. If equations are written into 
the manuscript in longhand, dubious-looking symbols 
should be identified with a penciled note. 


References should be listed on a separate sheet numbered 
in the order in which they occur in the text; or, if there are 
more than 10, in alphabetical order according to author. 
The listing should include author, title, and source (if the 
source is a magazine, the volume, number, month, year, 
and complete page numbers should be listed; if the source 
is a book, the publisher, place and date of publication, and 
the page numbers of the reference). If reference is made 


to a self-contained publication, the author, title, publisher, 
place of publication, and date should be given. Within 
the text, references should be indicated by Arabic numbers 
in brackets to correspond to the numbered list. 


Footnotes should be numbered consecutively in Arabic 
numerals and indicated in the text by superscripts. Each 
should be typed at the bottom of the page on which the 
footnote reference occurs. 


Illustrations.—A_ list of legends for the illustrations 
should be typed (double spaced ) on a separate sheet. Each 
illustration should be numbered in the margin or on the 
back outside the image area. To fit into the Review page, 
illustrations must take a reduction not to exceed 314 inches 
by 9 inches (column size) or 714 inches by 9 inches (page 
size). Map bases should show only political and conti- 
nental boundaries and latitude and longitude lines, unless 
data are to be plotted, when station circles will also be 
needed. Usually the less unnecessary detail in the back- 
ground the better will be the result from the standpoint of 
clear reproduction. Line drawings and graphs should 
also be uncluttered with fine background grids unless the 
graph demands very close reading. 


Photographs should be sharp and clear, with a glossy 
surface. Bear in mind that marks from paper clips or 
writing across the back will show up in the reproduction. 
Drawings and photographs should be protected with 
cardboard in mailing. 
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